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ABSTRACT 


Hamilton-type  geoacoustic  models  were  developed  for  Area 
Foxtrot,  a  NUWC  test  bed  for  emerging  active  sonar  systems  where 
the  surface  sediment  type  is  highly  spatially  variable. 

Reverberation  levels  (RL)  were  modeled  using  the  FEPE  propagation 
model  to  augment  the  GSM  propagation  model  because  the  bottom 
loss  model  in  GSM  did  not  estimate  transmission  loss  (TL) 
accurately  in  shallow  water.  FEPE  estimates  reveal  there  is  over  a 
15  dB  difference  between  TL  for  sand  and  that  for  silt-clay 
sediments  in  Area  Foxtrot.  The  comparison  between  modeled  RL  and 
measured  RL  (from  a  1991  ASW  exercise  conducted  by  NUWC)  enabled 
bottom  scattering  strength  kernels  to  be  developed  for  Area  Foxtrot. 
Bottom  scattering  strength  was  found  to  be  a  function  of  sediment 
type.  Hard  sand  sediment  has  a  bottom  scattering  strength  which 
obeys  Lambert’s  law  (sin2  (e))  while  that  of  silt-clay  sediment  is 

consistent  with  volume  scattering  (sin  (0)).  The  RLs  in  Area  Foxtrot 
are  azimuth-dependent  and  are  a  function  of  TL  and  bottom 
scattering  strength  (and  hence  bottom  sediment  type).  Sonar  beams 
steered  towards  the  hard  sand  show  higher  RLs  than  for  silt-clay. 
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I.  INTRODUCTION 


The  submarine  threat  to  Naval  forces  has  moved  from  deep 
water  environments  to  shallow  water.  Understanding  and  exploiting 
the  unique  acoustic  characteristics  of  coastal  waters  is  paramount 
to  our  national  security  and  the  success  of  Naval  forces  operating  in 
shallow  water. 

Shallow  water  environments  are  the  most  difficult 
environments  to  model  geoacoustically  because  of  the  variability  of 
the  bottom  sediment  type.  Not  only  are  the  surface  sediments 
laterally  inhomogeneous,  but  they  are  vertically  inhomogeneous  as 
well.  A  detailed  geologic  study  of  an  area  is  necessary  before  a  good 
geoacoustic  model  can  be  developed.  Geoacoustic  models  based  on 
sparse  geologic  data  are  of  little  value.  The  variability  of  the  ocean 
bottom  must  be  recognized  and  incorporated  into  the  Navy’s 
performance  prediction  models  to  accurately  estimate  transmission 
loss  and  reverberation  levels.  Failure  to  understand  the  variability 
of  ocean  bottom  sediments  means  inaccuracy  in  modeling  acoustic 
propagation  in  shallow  water  regimes.  This  in  turn  means 


1 


inaccuracy  in  modeling  sonar  performance  locally  and  in  potential 
conflict  areas. 

“The  littoral  (coastal)  region  is  frequently  characterized 
by  confined  and  congested  water  and  air  space  occupied 
by  friends,  adversaries,  and  neutrals--making 
identification  profoundly  difficult.  This  environment 
poses  varying  technical  and  tactical  challenges  to  naval 
forces.  It  is  an  area  where  our  adversaries  can 
concentrate  and  layer  their  defenses.  In  an  era  when 
arms  proliferation  means  some  Third  World  countries 
possess  sophisticated  weaponry,  there  is  a  wide  range  of 
potential  challenges.  For  example,  an  adversary’s 
submarines  operating  in  shallow  waters  poses  a 
particular  challenge  to  naval  forces. ..Mastery  of  the 
littoral  should  not  be  presumed.  It  does  not  derive 
directly  from  command  of  the  high  seas.  It  is  an 
objective  which  requires  our  focused  skills  and 
resources.”  - — excerpts  from  the  Navy  and  Marine  Corps 
White  Paper,  ...From  the  Sea  (O’Keefe  et  al.,1992). 

Many  potential  advisarial  nations  are  rapidly  acquiring 

/expanding  their  submarine  fleets.  Russia  has  agreed  to  sell  to  Iran 

three  (possibly  four)  Kilo-class  diesel  submarines,  the  first  of 

which  was  formally  commissioned  as  the  Taregh  (named  after  a 

famous  Muslim  warrior)  in  late  1992.  Recently  China  has  purchased 

a  Russian  Kilo-class  submarine.  Since  1985,  other  countries  which 

have  imported  Russian  Kilo-class  submarines  include  Algeria  (two). 
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Poland  (one),  Romania  (one),  and  India  (eight).  This  world-wide 
proliferation  of  stealthy  undersea  platforms  truly  supports  the  need 
for  continued  research  in  the  field  of  shallow  water  geoacoustics. 

Much  work  has  been  completed  in  geoacoustics  as  it  relates  to 
USW  (undersea  warfare)  in  deep  water.  In  contrast,  much  remains  to 
be  done  in  shallow  water  environments  especially  when  using  the 
active  sonar  systems  needed  to  prosecute  the  silent  diesel 
submarine.  Because  of  the  likelihood  of  multiple  bottom 
interactions  from  active  sonars,  operations  in  shallow  water  forces 
us  to  better  understand  the  influence  that  the  ocean  bottom  exerts 
on  acoustic  propagation.  If  the  geoacoustics  and  geophysics  of  the 
bottom  are  not  modeled  correctly,  each  bottom  interaction  will 
“amplify”  the  modeling  error  throughout  the  path  of  the  acoustic 
signal.  Even  small  errors  in  modeling  the  bottom  can  translate  into 
signal  excess  mis-calculations  of  many  decibels  in  magnitude. 

Better  detailed  geoacoustic  models  are  required  to  completely 
understand  an  acoustic  environment  and  to  quantify  and  qualify  the 
geophysics  involved  to  accurately  predict  the  performance  of 
current  and  emerging  sonars  in  shallow  water. 
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A.  REVERBERATION  DATA  BASE 


Reverberation  data  used  in  this  study  are  part  of  a  data  base 
collected  during  a  series  of  shallow  water  ASW  tests  conducted  3- 
13  September  1991  (Fisch  et  al.,  1991)  and  were  supplied  by  the 
Naval  Undersea  Warfare  Center  (NUWC),  New  London,  CT.  The  data 
supplied  to  the  Naval  Postgraduate  School  (NPS)  consisted  of  88 
pings  of  measured  reverberation  levels.  Half  of  the  pings  were 
continuous  wave  (CW),  500  millisecond  duration  waveforms  and  the 
other  half  were  coded  pulse  (CP),  500  millisecond  duration 
waveforms.  Each  ping  consisted  of  24  beams  with  a  total 

ft 

beamwidth  of  120  degrees.  The  depression  angle  of  the  sonar  was  3 
degrees.  This  study  only  used  the  CW  (500)  waveform  data. 

The  reverberation  data  were  collected  using  the  AN/SQS-53C 
sonar  onboard  USS  STUMP  (DD-978)  in  accordance  with  an 
established  test  plan  (NUSC,  1991).  The  data  used  in  this  study  were 
collected  on  7  September  1991  during  Run  18  of  the  ASW  tests.  Run 
18  was  conducted  in  64  m  of  water  with  a  negative  (downward) 
refracting  sound  speed  profile.  For  a  complete  environmental  report 
of  the  conditions  during  the  ASW  tests  see  Gallagher  (1991). 
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The  1991  ASW  tests  and  reverberation  data  collection  efforts 


were  conducted  in  Area  Foxtrot,  NUWC’s  test  bed  for  emerging 
sonars.  Area  Foxtrot  is  an  area  on  the  continental  shelf  of  the  U.S. 
eastern  seaboard  and  lies  approximately  50  nmi  south  of  Montauk 
Point,  Long  Island,  NY. 

B.  OBJECTIVES  OF  THIS  THESIS 
The  objectives  of  this  thesis  are: 

1.  Examine  the  current  geoacoustic  model  for  Area  Foxtrot  and,  if 
necessary,  develop  a  new  Hamilton-type  geoacoustic  model  for  the 
area. 

2.  Use  the  Generic  Sonar  Model  Multi-Path  Expansion  System 
propagation  model,  currently  used  by  NUWC  for  transmission  loss 
(TL)  and  reverberation  level  (RL)  modeling,  and  evaluate  the  model’s 
usefulness  in  a  shallow  water  environment. 

3.  Compare  measured  RL  and  modeled  RL  to  develop  a  bottom 
scattering  strength  kernel  for  Area  Foxtrot. 
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II.  THEORY  AND  BACKGROUND  OF  HAMILTON-TYPE 
GEOACOUSTIC  MODELS 


A  geoacoustic  model  of  the  seafloor  defines  the  thickness  and 
geophysical  properties  of  the  unconsolidated  sediment  and  rock 
strata  of  the  ocean  bottom.  Geoacoustic  models  are  used  in 
modeling  the  interaction  of  acoustic  energy  with  the  sea  bottom. 
Parameters  of  a  geoacoustic  model  consist  of  measured, 
extrapolated,  and  calculated  values  of  those  geophysical  properties 
important  in  estimating  the  propagation  of  acoustic  energy  (sound) 
when  it  interacts  with  the  seafloor.  Geoacoustic  models  guide 
theoretical  underwater  acoustic  studies,  help  reconcile  acoustic 
experiments  at  sea  with  theory,  and  aid  in  predicting  the  effects  of 
the  seafloor  on  acoustic  energy  (Hamilton,  1987).  The  general 
subject  of  geoacoustic  modeling  has  been  summarized  by  Hamilton 
(1980). 

Data  used  to  construct  a  geoacoustic  model  include 
bathymetric  charts,  echo-sounding  records  e.g.,  side-scan  sonars, 
analyses  from  seismic  reflection  and  refraction  records,  acoustic 
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reflection  records,  geologic  data  derived  from  sediment  or  rock 
grabs  and  cores,  and  water  mass  sound  speed  data.  A  wealth  of  such 
information  can  be  found  in  the  geologic  and  geophysical  literature 
and  through  today’s  Internet  resources. 

A  geoacoustic  model  should  contain  at  least  six  geoacoustic 
properties  for  each  sediment  or  rock  layer  to  be  of  value  to  a 
geophysicist  or  acoustician.  These  properties  are  1)  the  thickness 
of  the  layer  and  locations  of  any  reflectors  within  the  layer,  2)  the 
compressional  wave  speed  at  the  top  of  the  layer  and  the  gradient 
within  the  layer,  3)  the  compressional  wave  attenuation  at  the  top 
of  the  layer  and  the  gradient  within  the  layer,  4)  the  shear  wave 
speed  at  the  top  of  the  layer  and  the  gradient  within  the  layer,  5)  the 
shear  wave  attenuation  at  the  top  of  the  layer  and  the  gradient 
within  the  layer,  and  6)  the  density  at  the  top  of  the  layer  and  the 
gradient  within  the  layer.  Shear  wave  speed  and  shear  attenuation 
may  or  may  not  be  required  based  on  the  acoustic  model  used  and 
relative  magnitudes  of  compressional  and  shear  wave  speeds. 
Geoacoustic  properties  must  be  quantified  as  accurately  as  possible 
down  to  the  depth  of  the  acoustic  basement. 
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It  is  rare  that  all  six  properties  have  been  measured  for  any 
given  location.  Therefore,  the  missing  properties  must  be  estimated 
by  extrapolation  or  calculation  from  similar  sediments  in  the  same 
type  of  depositional  environment  but  often  at  distant  locations. 

Geoacoustic  modeling  in  shallow  water  environments  poses 
some  unique  challenges  and  difficult  problems.  Paramount  among 
these  are  the  vertical  and  horizontal  variability  of  sediment  and 
rock  types  and  the  consequent  variability  of  geoacoustic  properties. 
Surface  sediment  type  charts  based  on  a  few  bottom  samples  are 
insufficient  to  adequately  model  a  geoacoustic  environment  because 
the  uppermost  sediment  layer  may  be  relatively  thin  and  show 
lateral  pinch-outs  or  feathering  to  expose  underlying  strata.  Many 
seafloor  areas  were  originally  thought  to  be  spatially  homogeneous, 
and  thus,  acoustically  benign  based  on  some  wide-spread  surface 
sediment  data.  Seismic  reflection  profiling  with  sufficient 
resolution  would  enable  researchers  to  chart  the  lateral  extent  and 
thickness  of  reflectors.  Note  that  most,  if  not  all,  environments 
once  considered  acoustically  benign  are  actually  quite  spatially 
variable. 
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/ 

The  bottom  composition  in  shallow  water  environments  can  be 
highly  spatially  variable  within  one  general  local  and  throughout  the 
shallow  water  regime.  In  some  areas  the  bottom  sediment  may  be 
sand,  silt,  clay,  shell  fragments,  or  biogenic  ooze,  or  some 
combination  of  these.  In  other  areas  the  bottom  may  be  exposed 
rock.  This  rock  may  be  sedimentary,  igneous  (intrusive  or 
extrusive),  and  of  course,  metamorphic.  Some  rocky  seafloors  are  so 
rough  due  to  lava  flows,  large-scale  faults  or  small-scale  fractures, 
that  acoustic  energy  is  scattered,  reflected,  refracted,  and 
diffracted  to  such  a  great  extent  that  it  is  not  possible  to  construct 
a  useful  geoacoustic  model.  In  this  case,  the  bottom  reflectivity  or 
bottom  loss  at  the  frequency  of  interest  must  be  measured  in  situ. 

The  geoacoustic  models  developed  in  this  study  were 
constructed  using  a  method  put  forth  by  Hamilton  (1980)  and 
Hamilton  and  Bachman  (1982).  This  method  uses  a  series  of 
regression  equations  derived  from  measurements  of  sediment  and 
rock  geophysical  parameters  as  they  relate  to  the  geoacoustic 
properties  (sound  speed,  attenuation,  etc.).  The  geophysical 
parameters  Include  porosity,  grain  size,  density,  bulk  modulus,  and 


others.  Hamilton  (1987)  showed  that  the  geophysical  properties  of  a 
sediment  type  obtained  at  a  site  from  one  ocean  are  about  the  same 
as  those  for  a  sediment  in  another  ocean  given  that  the  two 
sediments  were  deposited  in  the  same  environment  and  under 
similar  depositional  conditions.  These  regression  equations  can 
therefore  be  used  world-wide  as  long  as  they  are  applied  to  the 
appropriate  geophysical  environment. 

Note  that  the  geoacoustician  must  have  a  good  understanding 
of  the  geologic  principles  governing  depositional  environments  and 
temporal  evolution  of  environments  in  order  to  temper  Hamilton's 
equations  and  apply  them  appropriately.  These  equations  should  not 
be  used  without  an  understanding  of  the  geological  sciences  and 
geologic  history  of  a  study  area.  For  instance,  the  Barents  Sea  is  a 
shallow  arctic  sea  environment.  But  direct  application  of 
Hamilton’s  equations  for  a  shallow  water  environment  to  the 
Barents  Sea  would  not  be  appropriate  considering  the  geologic 
history  of  the  area.  The  Barents  Sea  was  recently  (geologically 
speaking)  scoured  clean  of  its  bottom  unconsolidated  sediments  by 
glaciation.  The  enormous  mass  of  the  glaciers  produced  a  great  deal 


of  overburden  pressure  on  the  underlying  strata  increasing  the 
densities  of  the  strata  more  than  would  have  been  expected  in  the 
shallow  water  environment.  Geology  must  be  considered  in  the 
construction  of  geoacoustic  models. 

Acoustic  wave  speeds  (compressional  and  shear)  in  sediments 
are  functions  of  the  elastic  properties  of  the  two-phase  sediment 
mass  (pore  water  and  mineral  structure).  Hamilton  (1971,  1972) 
concluded  that  the  elastic  properties  of  water-saturated  sediments 
could  be  expressed  though  the  Hookean  elastic  equations  unless 
attenuation  is  considered.  In  this  latter  case,  linear  viscoelastic 
models  and  equations  are  recommended.  Sediment  porosity  and  grain 
size  effect  sound  speeds  only  through  their  effect  on  the  elastic 
properties  of  the  sediment. 

The  compressional  wave  speed  is  inversely  proportional  to 
porosity  and  is  directly  proportional  to  density.  The  compressional 
wave  speeds  in  a  low  porosity,  course-grained  sand  can  exceed  1800 
m/s  whereas  in  a  high  porosity  silt-clay  it  can  be  less  than  1500 
m/s  (sometimes  less  than  the  sound  speed  in  the  bottom  water) 
(Figure  1).  The  relationship  between  compressional  wave  speed  and 
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porosity  has  received  much  attention  because  porosity  is  an  easily- 
measured  sediment  property  and  the  relationship  usually  permits  a 
reasonable  prediction  of  sediment  acoustic  wave  speeds  (Hamilton, 
1987).  The  high  correlation  between  porosity  and  compressional 
wave  speed  occurs  because  it  is  the  volume  of  interstitial  pore 
water  that  largely  determines  the  compressional  wave  speed  in  the 
sediment.  This  is  a  result  of  the  relatively  high  compressibility  of 
water  compared  to  that  of  the  composite  mineral  grains  of  the 
sediment.  The  net  effect  of  the  varying  volumes  of  the  water  and 
mineral  grains  (and  hence  density)  results  in  many  continental 
terrace  silt-clays  and  most  abyssal  hill  and  plain  high  porosity 
sediments  having  compressional  wave  speeds  less  than  that  of  the 
overlying  bottom  water  (Hamilton,  1987).  Note  that  from  in  situ 
measurements,  compressional  wave  speeds  in  the  higher  porosity 
sediments  are  usually  no  more  than  three  percent  less  than  in  the 
bottom  sea  water  at  the  same  temperature  and  pressure. 

The  laboratory-derived  compressional  wave  speed  ratios 
between  the  sediment  and  water,  summarized  by  Hamilton  (1980) 
and  Hamilton  and  Bachman  (1982),  can  be  used  to  predict  the  in  situ 


sediment  surface  compressional  wave  speed  because  the  ratio  is  the 
same  in  the  laboratory  or  at  any  bottom  depth  in  the  ocean. 

However,  these  laboratory  values  must  be  corrected  for  depth  to  the 
seafloor.  To  make  the  correction,  the  laboratory  compressional 
wave  speed  ratio  is  multiplied  by  the  bottom  water  sound  speed. 

This  has  the  same  effect  as  making  full  temperature  and  pressure 
corrections  to  the  laboratory  measurements. 

The  relationship  between  sediment  density  and  compressional 
sound  speed  is  similar  to  that  for  porosity  because  of  the  linear 
relationship  between  density  and  porosity. 

The  sediment  mean  grain  size  from  all  types  of  bottom 
environments  and  the  percent  clay-sized  grains  from  deep  sea 
environments  were  determined  in  early  works  to  be  important 
indices  to  compressional  wave  speeds  (Hamilton,  1970;  Schreiber, 
1968;  Horn  et  al.,  1968;  Sutton  et  al.,  1957;  Hamilton  et  al.,  1956). 
Hamilton  and  Bachman  (1982)  verified  these  earlier  studies  and 
showed  that  the  acoustic  wave  speed  is  directly  proportional  to 
grain  size.  The  relationship  between  grain  size  and  acoustic  wave 
speed  (Figure  2)  is  especially  important  because  grain  size  analyses 
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can  be  made  from  dried  sediment  samples  for  which  density  (in 
situ),  porosity,  and  sound  speed  can  not  be  made.  Grain  size  can  also 
be  measured  from  simple  sediment  grabs  or  cores.  Additionally, 
much  data  exist  on  grain  size  in  the  geologic  literature  which  can  be 
used  as  indices  to  acoustic  properties. 

Grain  size  affects  the  acoustic  wave  speed  through  the 
properties  of  density  and  porosity.  However,  other  factors  such  as 
grain  shape,  mineralogy  and  the  presence  of  microcrystalline 
fractures  may  influence  the  structural  rigidity  and  therefore  the 
acoustic  wave  speed  of  the  sediment  in  minor  ways. 

The  calculation  of  compressional  wave  speed  can  be  performed 
using  the  porosity  (or  density)  or  grain  size  of  a  sediment.  Hamilton 
(1980)  and  Hamilton  and  Bachman  (1982)  have  compiled  a  series  of 
regression  equations  developed  from  in  situ  and  laboratory 
measurements  which  relate  these  properties  to  compressional  wave 
speed.  This  study  used  the  grain  size  regression  equations  to 
calculate  the  compressional  wave  speed  at  the  surface  of  the 
sediment  layer  since  these  data  were  most  available.  Verification 
of  the  calculated  compressional  sound  speeds  was  made  using 


available  sediment  density  data.  The  sediment  data  were  provided 
by  the  University  of  Rhode  Island  (Shoukimas,  1995). 

It  is  important  to  note  that  Hamilton’s  (1980)  and  Hamilton 
and  Bachman’s  (1982)  regression  equations  were  forced  to  fit 
through  the  acoustic  wave  speed  at  the  sediment  surface.  Acoustic 
wave  speed  equations  derived  from  sonobuoy  data  are  based  on  the 
mid-point  interval  velocity  from  the  first  layer  of  sediment.  No 
information  is  usually  available  from  sonobuoy  data  for  sound  speed 
at  or  near  the  surface  of  the  layer.  In  this  case,  the  intercepts  of 
regression  equations  for  the  top  of  the  sediment  layer  and  the  sound 
speed  gradient  from  the  surface  are  apt  to  be  statistical  artifacts 
(Hamilton,  1987).  A  specific  example  is  found  in  the  vast  amount  of 
published  equations  from  Lamont-Doherty  sources.  These  equations 
are  mostly  presented  as  linear  equations  and  the  intercepts  at  the 
seafloor  are  usually  statistical  artifacts.  The  equations  usually  can 
not  be  used  as  presented  to  calculate  instantaneous  layer  acoustic 
wave  speeds,  but  must  be  re-computed  with  realistic  seafloor 
acoustic  wave  speeds.  This  problem  is  documented  by  Lamont- 
Doherty  scientists  (Houtz,  1974;  Houtz  et  al.,  1968)  who  designed 


their  equations  for  primarily  computing  true  sediment  and  rock  layer 
thicknesses. 

Extrapolation  of  the  compressional  wave  speed  at  the  surface 
of  the  sediment  to  any  depth  within  the  layer  is  also  accomplished 
with  regression  equations  provided  by  Hamilton  (1979,  1985).  He 
used  laboratory  measurements  of  the  Shell  Development  Company  to 
compute  acoustic  wave  speed-depth  functions  in  sand  sediments. 
Hamilton  found  that  acoustic  wave  speed  gradients  in  sands  could  be 
very  high  relative  to  that  for  silt-clays.  For  example,  sands  could 
have  acoustic  wave  speed  gradients  twenty  times  that  for  silt-clays 
in  the  uppermost  1  -  3  m  of  the  sediment.  However,  the  acoustic 
wave  speed  gradient  for  silt-clay  is  variable  and  is  directly  related 
to  the  sedimentation  rate  and  diagenesis  of  the  sediment  (Hamilton, 
1987).  Diagenesis  refers  to  the  process  involving  physical  and 
chemical  changes  in  sediment  after  deposition  that  converts  it  to 
consolidated  rock.  For  silt-clays,  minimum  acoustic  wave  speed 
gradients  are  expected  in  and  near  areas  of  the  continental  terrace 
where  the  deposition  rate  is  rapid  (Hamilton  et  al.,  1974;  Houtz, 

1974;  Hamilton  et  al.,  1977;  Houtz,  1977;  Houtz,  1978;  Hamilton, 


1979;  Houtz,  1980;  Bachman  et  al.,  1983).  Here  again  geology  and 
geologic  history  of  a  study  area  is  required  to  accurately  estimate 
acoustic  wave  speed  gradients  in  sediments.  Note  that  the  greatest 
range  in  sediment  acoustic  wave  speeds  is  in  continental  terrace 
(and  hence  shallow  water)  environments. 

Once  the  compressional  wave  speed  has  been  determined 
throughout  the  depth  of  the  entire  sediment  layer,  the  compressional 
attenuation  for  the  layer  must  be  determined.  When  acoustic  energy 
is  refracted  through  a  sediment  layer  or  reflected  from  impedance 
contrasts  from  within  a  sediment  or  rock  layer,  some  energy  is 
returned  back  into  the  water  column.  The  amount  of  returned  energy 
is  a  function  of  the  distance  the  energy  has  traveled  through  the 
layer  and  the  attenuation  of  acoustic  energy  within  the  layer. 
Therefore,  the  compressional  wave  attenuation  at  the  surface  of  the 
layer  and  the  attenuation  gradient  within  the  layer  are  required 
parameters  of  a  geoacoustic  model. 

Hamilton  (1972,  1987)  determined  the  compressional  wave 
attenuation  in  sediments  as  a  function  of  acoustic  frequency.  It  is 
important  to  differentiate  between  intrinsic  attenuation  and 
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effective  attenuation.  Intrinsic  attenuation  is  the  loss  of  acoustic 


energy  through  conversion  to  heat  and  reduction  by  friction  between 
adjacent  sediment  grains.  Intrinsic  attenuation  is  the  type  usually 
measured  in  the  laboratory.  Effective  attenuation  includes  intrinsic 
attenuation  plus  attenuation  due  to  gas  bubbles,  animal  burrows, 
cobbles  and  large  stones,  and  large  shell  fragments  in  the  sediment.. 
Other  important  energy  loss  mechanisms  considered  part  of 
effective  attenuation  include  1)  losses  at  reflectors  and  refractors 
including  conversion  between  body  wave  types  (compressional  to 
shear)  and  conversion  of  body  waves  to  interface  or  surface  waves, 

2)  diffraction  and  scatter  due  to  reflector  roughness  and  curvature, 

3)  scattering  from  inhomogeneities  within  the  strata,  and  4) 
apparent  loss  due  to  inter-bed  multiples  which  add  destructively 
(Schoenberger  et  al.,  1974,  1975;  Scanlon  1985).  The  total  of  all 
energy  loss  mechanisms  (excluding  spreading  loss)  is  called 
‘effective  attenuation’.  It  is  effective  attenuation  which  the 
geoacoustician  must  use  to  fully  describe  the  energy  loss  within  a 
sediment  or  rock  layer.  Hamilton  (1972,  1987)  provides  measured 
attenuation  data  for  sediments  from  various  depositional 


environments.  Note  that  Hamilton’s  measured  data  include  effective 


attenuation  (of  which  intrinsic  attenuation  is  a  part). 

Sediment  type  and  its  physical  properties  greatly  effect  the 
degree  of  attenuation  in  a  layer.  Sands  can  attenuate  acoustic 
energy  more  than  ten  times  that  of  silt-clays  (Figure  3).  It  is 
therefore  difficult  to  predict  attenuation  in  a  sediment  layer 
without  a  detailed  knowledge  of  the  sediment  composition  at  the 
surface  and  with  depth.  Although  the  compressional  attenuation  for 
sand  sediments  is  greater  than  that  for  silt-clay  sediments,  per  unit 
distance  traveled,  the  overall  total  attenuation  for  a  sand  layer  is 
generally  much  less  than  that  for  silt-clay  layers  because  less  total 
distance  is  traveled  by  acoustic  energy  in  sand.  This  is  due  to  the 
greater  compressional  wave  speed  gradient  for  sand  than  for  silt- 
clay.  Sand  is  therefore  a  more  efficient  refractor  of  acoustic  energy 
and  returns  acoustic  energy  back  to  the  water  column  with  less 
internal  travel  distance  than  for  silt-clay. 

A  compilation  of  attenuation  data  by  Hamilton  (1972)  and 
Hamilton  and  Bachman  (1982)  suggest  that  attenuation  is  related  to 
frequency  through  a  constant  as  in  the  equation 
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ttp  =  kp  f 


(1) 


where  ap  is  the  compressional  wave  attenuation  in  dB/m 
kp  is  an  attenuation  constant  (dB/m-kHz) 
f  is  the  acoustic  frequency  (kHz) 

If  attenuation  is  dependent  on  the  first  power  of  frequency  as  in 
Equation  1,  the  constant  kp  must  vary  with  sediment  or  rock  type  and 
be  a  function  of  porosity  and  grain  size.  Given  either  porosity  or 
grain  size,  an  estimate  of  the  compressional  wave  attenuation  can 
be  found  using  the  figures  and  equations  provided  by  Hamilton 
(1972).  Hamilton  (1987)  gives  compressional  wave  attenuation 
estimates  for  .finer-grained,  higher  porosity  sediments  than  those  in 
Hamilton  (1972). 

Extrapolating  the  compressional  wave  attenuation  from  the 
sediment  surface  to  deeper  depths  in  the  sediment  layer  is  perhaps 
the  most  difficult  aspect  of  geoacoustic  modeling.  Silt-clay 
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sediments  (fine-grained)  show  an  increase  in  attenuation  with  depth 
down  to  about  500  m  after  which  the  attenuation  decreases  with 
increased  depth  (Figure  3).  Attenuation  values  for  silt-clays  can 
vary  with  depth  from  very  low  attenuation  values  (Mitchell  and 
Focke,  1980;  Frisk  et  al.,  1981,  1986)  to  much  higher  values 
(Neprochnov,  1971;  Jacobson  et  al.,  1981,  1984)  which  are 
associated  with  clay-sands  and  shales.  It  is  likely  that  lithification 
of  the  silt-clays  into  shale  is  responsible  for  the  increased 
attenuation  values. 

The  differences  in  attenuation  with  depth  from  the  sources 
above  are  probably  due  to  differences  in  sediment  composition, 
consolidation  and  lithification,  diagenesis,  and  stratigraphy.  These 
are  many  of  the  same  factors  which  cause  variations  in  sound  speed. 

For  sands,  the  attenuation  changes  with  depth  are  a  bit  more 
predictable.  Sands  show  a  decrease  in  attenuation  with  increasing 
depth  at  all  depths.  Hamilton  (1976,  1987)  suggested  that 
attenuation  in  sands  decreases  with  about  -1/6  the  power  of  the 
overburden  pressure.  He  further  suggested  that  to  predict  the  value 
of  attenuation  in  sand,  one  selects  a  value  of  kp  (see  Equation  1)  at 
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the  surface  and  then  decrease  this  value  at  a  rate  of  (depth)-i/6.  a 
constant  attenuation  with  depth  certainly  can  not  be  assumed 
(Mitchell  and  Focke,  1983). 

As  with  sound  speeds,  attenuation  values  in  each  sediment  and 
rock  layer  is  apt  to  be  unique.  The  best  that  can  be  done  for 
prediction  of  effective  attenuation  and  its  gradient  is  to  gather  data 
from  similar  sediment  types  and  geomorphic  areas,  examine  the  data 
statistically,  and  relate  this  Information  between  study  areas.  No 
theory  or  theoretical  approach  (e.g.,  Biot-Stoll)  can  be  used  to 
accurately  compute  the  effective  attenuation  with  depth  because  of 
the  variability  of  energy  loss  mechanisms  beneath  the  surface. 

Currently,  the  geoacoustician  can  begin  with  an  estimate  of 
attenuation  at  the  sediment  surface  and  then  vary  the  attenuation 
value  with  depth  according  to  an  estimate  of  the  amount  of  sand  and 
silt-clay  in  the  geologic  section.  For  guidance,  the  data  of 
Hamilton  (1976)  can  be  followed  if  there  are  many  sand-silt  layers. 

If  much  less  (or  no)  amounts  of  sand-silt  are  expected,  the  data  of 
Mitchell  and  Focke  (1980)  can  be  used. 
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Hamilton  (1987)  agrees  with  Stoll  (1983,  1985)  that  it  is  not 
possible  to  compute,  with  any  theory,  the  effective  attenuation  as  a 
function  of  depth  in  thick  strata.  Theories  can  relate  frequency  and 
burial  depth  of  sediments  to  intrinsic  attenuation  but  not  to 
effective  attenuation. 
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III.  PRELIMINARY  GEOACOUSTIC  MODEL  DEVELOPED  BY  NPS 
AND  ENVIRONMENTAL  SUMMARY 


NUWC’s  test  bed  for  emerging  sonar  systems  is  Area  Foxtrot. 
NUWC  considers  Area  Foxtrot  “significant  in  that  it  has  been  shown 
to  have  bottom  composition  and  seasonal  sound  speed 
characteristics  similar  to  several  potential  regional  conflict  areas” 
(Christian  and  Capizzano,  1994).  For  this  reason  the  historically 
accepted  geoacoustic  model  for  Area  Foxtrot  (Cole  et  al.,  1994)  was 
examined  by  NPS  to  seek  areas  of  possible  model  improvement.  This 
current  geoacoustic  model  for  Area  Foxtrot  (discussed  below)  has 
been  used  for  both  past  and  recent  acoustic  studies  An  improved 
geoacoustic  model  developed  in  this  study  is  based  on  Hamilton’s 
(1970,  1971,  1972,  1976,  1979,  1980,  1985,  1987)  and  Hamilton  and 
Bachman’s  (1982)  geoacoustic  modeling  techniques  that  are  now 
well-established.  These  models  will  incorporate  Hamilton’s 
geoacoustic  parameters  derived  from  a  variety  of  geoacoustic  and 
geologic  data. 
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A.  BACKGROUND  OF  NUWC’S  CURRENT  GEOACOUSTIC  / 
BOTTOM  LOSS  MODEL 

In  May  1965  TL  measurements  at  3.5  kHz  were  made  in  Area 
Foxtrot  by  the  Naval  Underwater  Sound  Laboratory  (Fisch  et  al., 

1992).  Analyses  at  that  time  indicated  that  the  measured  bottom 
loss  was  1.3  dB  per  bounce  at  a  grazing  angle  of  11  degrees.  An 
early  ray-tracing  program  was  used  to  estimate  the  spreading  loss, 
and  with  questionable  assumptions  made  about  volume  attenuation,  a 
bottom  loss  of  1.3  dB  per  bounce  was  determined.  This  value  was 
found  to  be  in  good  agreement  with  a  bottom  loss  value  extracted 
from  a  geoacoustic  model  which  was  based  on  sediment  core 
information  from  the  TL  measurement  site  (Cole  et  al.,  1994).  The 
calculated  bottom  loss  curve  (Cole  et  al.,  1994)  is  based  on  a  half¬ 
space  mathematical  model  which  includes  attenuation.  It  used  a 
sediment  density  (from  core  data)  of  2.0  g/cm3  (Cole,  1956)  and  a 
depth-independent  compressional  wave  speed  of  1554  m/s  using  the 
dispersion  method  (Barakos,  1962).  These  values  are  in  stark 
disagreement  with  the  vast  number  of  geoacoustic  measurements 
made  by  Hamilton  (1976,  1979,  1980).  The  bottom  attenuation. 
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inferred  from  the  above  measurements,  was  0.45  dB/ft.  Assuming  a 
linear  dependence  of  attenuation  with  frequency,  this  corresponds  to 
an  attenuation  parameter  (kp)  of  0.422  dB/m-kHz.  Table  1  lists  the 
geoacoustic  parameters  used  in  the  model.  The  sediment 
compressional  wave  speed  based  on  the  above  values  (Cole  et  al., 
1994)  is  plotted  in  Figure  4  which  also  incorporates  the  water  sound 
speed  profile  measured  for  Run  18.  The  calculated  bottom  loss  curve 
is  frequency  independent  since  the  half-space  geoacoustic  model 
used  a  constant  attenuation  coefficient. 

A  bottom  loss  curve,  calculated  from  the  geoacoustic  model 
(see  Figure  1  of  Cole  et  al.,  1994),  was  used  as  an  input  to  the  GSM 
model  for  active  sonar  performance  predictions  by  NUWC 
researchers  (Cole  et  al.,  1994).  That  study  assumed  a  Lambert’s  law 
of  diffuse  backscattering  at  the  ocean  bottom.  The  value  for  the 
constant  “k”  used  to  determine  the  bottom  backscatter  strength  at 
3.5  kHz  was  derived  from  a  1965  Area  Foxtrot  measurement  study  to 
be  -19.4  dB.  The  k  variable  represents  the  quantity  10  log  p  in  the 
Lambert’s  Law  relationship. 
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Sb  =  10  log  p.  +  10  log  sin2  e 


(2) 


where  Sb  is  bottom  backscattering  strength 

|x  is  a  site-specific  empirical  value  and  is  a  function  of 
bottom  composition  and  roughness,  and 
0  is  the  incident  angle. 

Cole  et  al.  noted  that  the  k  value  for  Area  Foxtrot  was  about  6 
dB  greater  than  the  empirical  value  derived  by  Mackenzie  (Cole  and 
Podeszwa,  1974;  Mackenzie,  1961).  The  value  of  k  measured  in  Area 
Foxtrot  is  reported  to  vary  between  -18  dB  and  -21  dB  (Bishop, 

1994).  This  is  in  contrast  to  Mackenzie’s  world-wide  average  k 
value  of  -29  dB.  A  brief  examination  of  Lambert’s  law  will  be  found 
later  in  this  text. 

Values  of  k  ranging  from  -18  to  -21  dB  is  a  good  initial 
estimate,  but  may  not  be  accurate  considering  the  geoacoustic  model 
developed  by  NFS  as  part  of  this  thesis.  The  current  backscattering 
strength  (k)  is  based  on  essentially  the  same  models  and  procedures 
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used  to  calculate  the  bottom  loss  curve  and  uses  the  assumption 
that  reverberation  decay  depends  directly  on  bottom  loss  (Cole  et 
al.,  1994).  Inaccurate  geoacoustic  parameters  in  the  current 
geoacoustic  model  (Table  1)  are  expected  to  degrade  the  accuracy  of 
the  calculated  bottom  loss  curve,  reverberation  decay  rate  and 
estimated  value  of  k.  Because  of  the  influence  of  the  geoacoustic 
model  on  these  parameters,  this  study  will  examine  the  geoacoustic 
parameters  of  part  of  Area  Foxtrot  to  validate  the  current 
geoacoustic  model  in  use  at  NUWC.  Since  the  reported  model 
estimates  of  TL  are  consistently  greater  than  TL  measurements  in 
Area  Foxtrot  (Cole  et  al.,  1994),  the  potential  for  an  error  in  the 
geoacoustic  model  is  considered  likely  since  sediment  parameters  in 
the  model  are  inconsistent  with  Hamilton’s  geoacoustic  models.  It 
will  also  be  shown  that  the  spatial  variability  of  the  surface 
sediment  type  in  Area  Foxtrot  is  highly  variable  and  a  single 
geoacoustic  model  applied  to  the  entire  area  is  not  adequate  for 
performance  prediction.  This  leads  to  the  conclusion  that 
geoacoustic  modeling  must  be  based  on  Hamilton’s  measurements 
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and  that  the  geoacoustic  model  must  be  spatially  varied  as  surface 
sediment  type  varies. 

B.  RECOMMENDED  CHANGES  TO  NUWC’S  CURRENT 
GEOACOUSTIC  MODEL 

The  current  geoacoustic  model  uses  a  surface  sediment 
compressional  wave  speed  of  1554  m/s  which  was  measured  in  Area 
Foxtrot  from  data  collected  in  1958  and  1959  (Barakos,  1962).  This 
compressional  wave  speed  is  much  too  slow  for  a  sand  sediment 
with  a  density  of  2.0  g/cm3  (Hamilton,  1980)  and  an  overall  mean 
grain  size  (i.e.,  spatially  averaged  from  the  mean  grain  size  of  each 
sediment  sample)  of  2.01  phi  units  (Brocher  and  Ewing,  1986)  in  a 
continental  terrace  depositional  environment  (Hamilton,  1982).  Note 
that  phi  size  is  -loga  of  the  grain  diameter  in  millimeters.  A 
compressional  wave  speed  of  1554  m/s  is  more  consistent  with  a 
clayey-silt  sediment  on  the  continental  terrace.  The  associated 
grain  phi  size  for  a  clayey-silt  is  7.13  phi  units.  An  obvious 
discrepancy  exists  since  the  overall  mean  phi  size  was  measured  as 
2.01  phi  units  (Brocher  and  Ewing,  1986)  and  the  sediment  was 
reported  to  be  sand  (Brocher  and  Ewing,  1986;  Cole,  1956)  and  not 
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clayey-silt.  The  attenuation  value  of  0.422  dB/m/kHz  as  measured 
by  NUWC  is  also  applicable  to  a  sand  sediment  and  is  much  too  high 
for  a  sediment  with  a  compressional  wave  speed  of  1554  m/s 
(Hamilton,  1980).  Based  on  these  considerations,  the  sediment 
compressional  wave  speed  is  likely  much  faster  than  1554  m/s. 
Hamilton’s  geoacoustic  measurements,  made  in  many  depositional 
environments,  never  show  sediments  to  have  such  a  high  density  and 
at  the  same  time  be  associated  with  the  low  compressional  wave 
speed  of  1554  m/s. 

The  dispersion  method  Barakos  used  to  determine  the  bottom 
sediment  compressional  wave  speed  may  not  in  fact  be  a  viable 
method  for  computing  the  speed  of  compressional  waves  (body 
waves).  Barakos  (1962)  used  a  stationary  tower,  Texas  Tower  #4, 
near  Area  Foxtrot  as  a  stationary  receiving  station  to  which 
information  was  relayed  via  cable  from  hydrophones  in  the  ocean. 
Ships  were  deployed  from  Texas  Tower  #4  in  five  radial  arms 
outward  and  away  from  the  tower  to  a  range  of  150  miles  (Figure  1 
of  Barakos,  1962).  One  pound  explosive  charges  were  dropped  by  the 
ships  at  various  ranges  as  they  headed  away  from  the  tower. 
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Pressure  waves  generated  by  the  explosions  were  received  by  the 
hydrophones  and  data  relayed  to  the  tower.  Dispersion 
characteristics  of  the  pressure  waves  were  studied  and  used  to 
calculate  the  compressional  wave  speed  in  the  bottom  sediments. 
The  average  compressional  wave  speed  along  all  five  radial  arms, 
each  of  which  extends  150  miles  in  different  directions,  was 
reported  to  be  1554  m/s.  However,  it  is  not  geologically  feasible  to 
use  this  method  over  the  large  surface  area  in  which  Barakos 
employed  the  technique,  about  45,000  square  miles,  and  then  take 
the  average  of  the  compressional  wave  speed  of  all  the  surficial 
sediments  and  apply  this  average  compressional  wave  speed  to  the 
entire  study  area.  The  study  area  includes  continental  terrace, 
continental  slope  and  submarine  canyon  geologic  settings,  each  with 
unique  sedimentological  processes  and  depositional  characteristics. 
Within  any  of  the  above  geologic  settings,  the  sediment 
compressional  wave  speed  will  vary  significantly.  Averaging  of 
sediment  compressional  wave  speed  between  and  throughout 
depositional  environments  yields  no  useful  compressional  wave 
speed  for  any  of  the  environments. 
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Perhaps  of  more  importance  is  the  fact  that  Barakos 
calculated  the  sediment  sound  speed  based  on  dispersion  of  the 
ground  wave.  Some  surface  waves  such  as  Rayleigh  and  Airy  waves 
undergo  dispersion.  Love  waves,  another  type  of  surface  wave,  do 
not.  Guided  waves  undergo  dispersion  as  well.  In  dispersive  waves, 
particle  motion  is  described  as  an  elliptical,  retrograde  orbit. 
Therefore,  the  speed  of  the  wave  depends  on  the  elastic  constants  of 
material  across  an  interface  and  not  merely  on  elastic  constant 
values  at  a  given  point  as  for  shear  and  compressional  body  waves. 
Barakos’  calculated  speeds  for  the  sediment  are  for  waves  at  a  given 
frequency  and  are  not  representative  of  compressional  wave  speeds 
within  the  sediment.  Compressional  and  shear  waves  show  such 
small  dispersion  that  special  investigations  to  determine  body  wave 
dispersion  have  failed  (Richter,  1958). 

C.  CONSTRUCTION  OF  A  HAMILTON  GEOACOUSTIC  MODEL  FOR 
AREA  FOXTROT 

Brocher  and  Ewing  (1986)  report  that  data  were  collected  by 
the  United  States  Geological  Survey  in  Area  Foxtrot  in  1982,  during 
which  a  12-channel,  high-resolution  seismic  reflection  study  was 
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performed.  These  data  were  collected  in  approximately  64  m  of 
water,  close  to  an  existing  single  channel  reflection  profile.  They 
found  the  seafloor  to  be  underlain  by  about  2000  to  3000  m  of 
sediment,  the  uppermost  of  which  consists  of  a  thin  (approximately 
30  m)  layer  of  unconsolidated  Pleistocene  and  younger  sediment. 

This  thin  sediment  layer  unconformably  overlies  a  middle  Tertiary 
deltaic  sequence  which  thickens  toward  the  continental  shelf  edge. 
The  30  m  surface  sediment  layer  consists  of  reworked  sediments 
from  older,  underlying  sediments  as  well  as  riverine  clastic  debris 
deposited  during  times  of  lower  sea  level.  Nine  sediment  samples 
were  collected  within  50  km  of  the  seismic  study  area.  Sediments 
recovered  were  coarse-to-fine  sand.  The  sands  coarsened  toward 
Long  Island;  the  overall  mean  grain  size  was  2.01  phi  units.  The 
density  of  the  sediment  was  found  from  a  core  sample  to  be  2.0 
g/cm3  (Cole,  1956)  and  is  consistent  with  the  density  of  a  fine-to- 
medium  grain  sand.  The  exact  location  of  this  core  sample  is  not 
currently  known. 

In  a  continental  terrace  depositional  environment  sediment 
grain  size  will  decrease  seaward  towards  the  continental  shelf;  this 


34 


is  consistent  with  the  sediment  samples  from  around  Area  Foxtrot 
(Brocher  and  Ewing,  1986).  Since  Run  18  (Fisch  et  al.,  1992),  the 
subject  of  our  backscattering  analysis,  took  place  in  the 
southeastern  part  of  Area  Foxtrot  and  nearest  to  the  continental 
slope,  the  grain  phi  size  should  be  greater  than  the  overall  mean  phi 
size  of  2.01  phi  units  reported  by  Brocher  and  Ewing  (1986)  (greater 
phi  size  means  smaller  grain  size  diameter).  The  phi  size  of 
sediment  in  the  Run  18  area  was  originally  thought  to  be  about  2.4 
phi  units  or  that  which  corresponds  to  a  fine  sand.  This  best-guess 
for  grain  size  was  made  with  little  geologic  data  and  no  exact  core 
or  sediment  sample  locations. 

However,  an  in-depth  data  search  for  sediment  grab  and  core 
data  has  revealed  much  geologic  data  exists  for  the  sediments  in 
Area  Foxtrot.  A  compilation  of  geologic  data  sets  for  Area  Foxtrot 
and  the  surrounding  area  was  obtained  through  a  data  search 
conducted  at  the  University  of  Rhode  Island  (Shoukimas,  1995).  The 
locations  of  where  sediment  samples  were  retrieved  within  Area 
Foxtrot  are  seen  in  Figure  5.  The  method  of  calculation  of 
compressional  wave  speed  is  discussed  below. 
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The  geoacoustic  models  developed  in  this  study  were 
constructed  using  Hamilton’s  (1972,  1976,  1979,  1982)  series  of 
equations.  One  of  Hamilton’s  (1982)  regression  equations  employed 
here  uses  the  sediment  mean  phi  size  to  calculate  the  compressional 
wave  speed  for  continental  terrace  unconsolidated  sediments: 

Vp  =  1952.5  -  86.26  {^)  +  4.14  (({>2)  (3) 

where  ^  is  mean  grain  size  in  phi  units 

Vp  is  compressional  wave  speed  in  m/s 

Equation  3  was  used  to  calculate  the  compressional  wave 
speeds  in  this  study.  Figure  6  is  a  contour  plot  of  the  mean  grain 
size  in  phi  units  in  and  around  Area  Foxtrot.  An  expected  general 
coarsening  of  grain  size  toward  the  continental  coast  (northwest)  is 
observed  with  smaller  grain  sizes  (silt  and  clay)  seen  on  the 
continental  slope  (southeast).  Figure  7  is  a  contour  plot  of  the 
compressional  wave  speeds  calculated  from  the  mean  grain  size  in 
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and  around  Area  Foxtrot.  The  spatial  variability  in  sediment  sound 
speed  mirrors  that  of  the  mean  grain  size.  From  Equation  3  it  can  be 
seen  that  relatively  small  changes  in  grain  size  can  translate  into 
significant  changes  in  compressional  wave  speed. 

Figures  3  and  4  show  the  variability  of  surface  sediment  type 
in  Area  Foxtrot.  Some  bottom  areas  are  covered  by  coarse  sand 
while  others  are  covered  by  silt-clay.  These  are  only  the  two 
extremes  of  surface  sediment  type.  There  are  areas  of  the  bottom 
that  are  covered  with  any  combination  of  sand  and  silt-clay. 
Therefore,  there  is  no  one  geoacoustic  model  which  correctly 
describes  all  of  Area  Foxtrot  as  originally  presumed  by  NUWC. 

Figure  8  depicts  the  tracks  and  approximate  location  of  the 
surface  ship  and  target  for  all  runs  of  the  September  1991  study 
(Cole  et  al.,  1994)  plotted  with  the  sediment  compressional  wave 
speed  in  shades  of  grey.  Although  the  contour  interval  is  fine  (50 
m/s),  It  can  be  seen  that  the  surface  ship  and  target  were  not  in  the 
same  acoustic  environment.  This  is  most  evident  for  Runs  1-8.  Runs 
11-18  show  the  least  variability  in  acoustic  environment  between 
the  location  of  the  ship  and  target. 
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The  degree  to  which  the  difference  in  acoustic  environment 
affects  the  transmission  loss  and  reverberation  level  between  the 
source  and  target  is  in  part  a  function  of  frequency  and  will  be 
addressed  in  future  research  efforts. 

Figure  9  shows  some  sediment  parameters  around  the  surface 
ship  track  of  Run  18.  Notice  the  variability  of  the  surface  sediment 
even  in  this  relatively  small  area.  Interpolating  the  data,  the 
sediment  directly  under  Run  18  is  likely  100%  coarse  sand  with  a 
first  approximation  of  compressional  wave  speed  of  about  1850  m/s. 
This  value  is  consistent  with  laboratory  values  for  a  coarse  sand  on 
the  continental  terrace  (Hamilton,  1982)  though  it  is  slightly  above 
Hamilton’s  mean  compressional  wave  speed  of  1836  m/s.  The 
standard  error  of  the  mean  of  the  laboratory  measurements  was  not 
determined. 

The  associated  average  in  situ  sediment/water  compressional 
wave  speed  ratio  for  a  coarse  sand  bottom  is  1.201  (Hamilton  and 
Bachman,  1982).  During  Run  18  the  water  sound  speed  immediately 
above  the  seafloor  was  1484  m/s,  measured  within  1.5  hours  of 
collecting  the  reverberation  data.  Based  upon  the  calculated 
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sediment  compressional  wave  speed  of  1850  m/s,  a  sediment  /water 
sound  speed  ratio  of  1.246  was  obtained.  This  ratio  is  slightly 
higher  than  the  above  value  determined  by  Hamilton  and  Bachman 
(1982)  for  a  continental  terrace  environment,  and  is  significantly 
different  from  the  sediment/water  sound  speed  ratio  of  1.047 
obtained  when  a  sediment  compressional  sound  speed  of  1554  m/s  is 
used  from  the  current  geoacoustic  model.  This  latter  ratio  is  far  too 
low  and  out  of  range  for  a  sandy  sediment.  Hence,  the  compressional 
wave  speed  of  1554  m/s,  used  in  previous  analyses  for  this  area,  is 
inconsistent  with  the  measured  geoacoustic  parameters  in  Area 
Foxtrot. 

The  relationship  between  grain  size  and  compressional  wave 
speed  is  especially  important  because  grain  size  analyses  can  be 
made  for  dried  sediment  for  which  measurements  of  density, 
porosity,  and  velocity  can  not  be  made.  Additionally,  much  data  on 
grain  sizes  exists  in  the  geologic  literature  which  can  be  readily 
used  as  Indices  of  acoustic  properties.  Hamilton’s  regression 
equation  for  compressional  wave  speed  as  a  function  of  grain  size  is 
as  good  or  slightly  better  than  the  regression  equation  based  on 
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sediment  density.  The  standard  error  of  estimate  for  mean  grain 
size  versus  sound  speed  in  continental  terrace  sediments  is  29  m/s 
(Hamilton  and  Bachman,  1982).  Hence,  the  compressional  wave 
speed  in  the  vicinity  of  Run  18  could  range  from  1821  m/s  to  1879 
m/s  vice  1554  m/s  as  used  in  the  current  geoacoustic  model. 

As  an  additional  validity  check  on  our  calculated  sediment 
compressional  wave  speed,  the  average  sediment/water  sound  speed 
ratio  of  1.201  (Table  2  of  Hamilton  and  Bachman,  1982)  can  be 
multiplied  by  the  sound  speed  of  the  bottom  water  in  our  study  area 
(1484  m/s)  which  yields  an  average  in  situ  sediment  compressional 
wave  speed  1782  m/s.  Multiplying  the  sediment/water  sound  speed 
ratio  by  sediment  compressional  wave  speed  has  the  effect  of 
correcting  laboratory  values  for  full  temperature  and  pressure 
corrections.  Given  this  rather  low  in  situ  speed  (1782  m/s)  and  the 
low  end  of  the  laboratory  speed  calculation  (1821  m/s),  it  is  likely 
that  the  true  compressional  sound  speed  lies  somewhere  between 
the  two.  However,  it  would  not  be  a  straight  average  of  the  two 
speeds.  Since  the  lithology  of  the  sediment  grab  nearest  to  Run  18  is 
a  coarse,  well-sorted  sand,  the  true  compressional  wave  speed  of 
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the  sediment  is  likely  much  closer  to  1821  m/s.  For  modeling 
purposes,  this  study  will  use  a  compressional  wave  speed  1817  m/s, 
similar  to  the  compressional  wave  speed  measured  in  situ  by 
Hamilton  (1972)  for  a  coarse  sand  in  34  m  of  water  off  the  coast  of 
San  Diego.  This,  again,  gives  another  check  on  the  validity  of  the 
calculations  for  sediment  compressional  wave  speed.  The 
compressional  wave  speed  given  above  is  in  relatively  good 
agreement  with  the  sediment  compressional  wave  speed  of  1750 
m/s  measured  during  the  seismic  survey  of  Area  Foxtrot  (Brocher 
and  Ewing,  1986).  It  must  be  noted  that  the  seismic  method  used  in 
the  1986  study  lends  itself  to  an  error  of  several  hundred  meters 
per  second  in  determining  sound  speed  in  the  sediment.  Table  2 
gives  the  elements  of  the  NFS  geoacoustic  model. 

To  construct  a  more  complete  geoacoustic  model,  changes  in 
the  sediment  compressional  wave  speed  (m/s)  as  a  function  of  depth 
in  the  sediment  column  due  to  overburden  pressure  were  calculated 
using  Hamilton’s  (1979)  regression  equation: 


where  K  is  a  constant  (s-i)  and 

D  is  depth  in  the  sediment  column  (m). 

To  calculate  K  for  our  area  of  study  consisting  of  a  coarse  sand 
bottom  the  following  expression  (Hamilton,  1979)  is  used: 

K  =  Vpo  /  0.05  0.015  (5) 

where  Vp©  is  the  sediment  compressional  wave  speed  at  the  surface 
and  the  value  0.05  brings  the  surface  compressional  wave  speed  to  a 
depth  of  5  cm  in  the  sediment  and  helps  to  correct  for  laboratory 
measurements.  If  speeds  calculated  in  the  laboratory  only  are 
considered,  K  would  be  1901.5  s-i  and  Vp  at  the  bottom  of  the 
Pleistocene  sediment  layer  (30  m)  would  be  2001  m/s.  Because  the 
determination  of  sediment  compressional  wave  speed  is  in  part 
based  on  Hamilton’s  (1972)  in  situ  measurements  of  speed  obtained 
at  depths  in  the  sediment  between  30  cm  and  60  cm,  a  K  value  of 
1817  s-1  is  deduced  and  the  initial  correction  for  depth  is  neglected. 
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The  final  equation  which  relates  sound  speed  to  depth  is  then 

Vp  =  1817  (D  0.015)  (6) 

and  the  resulting  compressional  wave  speed  at  the  bottom  of  the 
sand  layer  (30  m)  is  1945  m/s.  As  seen  in  Equations  5  and  6,  the 
increase  of  speed  with  depth  is  not  linear. 

Attenuation  at  3.5  kHz  was  determined  using  Hamilton’s 
(1972)  in  situ  attenuation  measurements  for  coarse  sand.  He 
documented  measurements  of  compressional  wave  speed  and 
attenuation  made  by  himself  as  well  as  from  others  found  in  the 
literature.  Of  note  are  three  important  relationships  between 
frequency  and  other  physical  parameters  which  Hamilton  observed 
from  his  analyses:  (1)  compressional  wave  attenuation  in  marine 
sediments  is  approximately  related  to  the  first  power  of  frequency 
in  sands,  muds,  and  sedimentary  rocks,  (2)  attenuation  can  be 
estimated  from  sediment  grain  size  (porosity)  data,  and  (3) 
attenuation  in  marine  sediments  typically  increases  with  depth  due 
to  the  reduction  of  sediment  porosity  to  a  point  where  overburden 
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pressure  becomes  the  dominant  effect.  Beyond  this  point 
attenuation  decreases  smoothly  as  a  function  of  depth  and  increased 
overburden  pressure  (Hamilton,  1976).  Hamilton  (1976)  incorporated 
other  attenuation  measurements  into  his  study  and  substantiated 
observation  (1)  above. 

At  3.5  kHz,  the  attenuation  for  a  coarse  sand  with  a 
compressional  wave  speed  of  1817  m/s  in  32  m  of  water  off  the 
coast  of  San  Diego  was  measured  in  situ  by  Hamilton  to  be  0.53 
dB/m-kHz  (kp  in  Hamilton’s  notation).  This  corresponds  to  an 

attenuation  coefficient,  ap,  of  1.76  dB/m  at  3.5  kHz.  Since  the 
coarse  sand  in  Area  Foxtrot  where  Run  18  took  place  is  in  deeper 
water  (64  m),  the  attenuation  may  be  slightly  less  due  to  a  small 
decrease  in  sediment  porosity.  But  because  the  sand  is  well-sorted 
and  there  Is  only  an  additional  30  m  more  water  overlying  the  Area 
Foxtrot  sediments  than  the  sediments  off  San  Diego,  the  reduction  in 
porosity  will  be  very  small.  Attenuation  may  be  closer  to  0.5 
dB/m-kHz.  This  attenuation  value  is  higher  than  that  of  0.422 
dB/m-kHz  (1.477  dB/m  at  3.5  kHz)  used  in  the  current  geoacoustic 
model. 
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Caution  must  be  taken  regarding  attenuation  coefficients. 
Hamilton  (1976)  defines  attenuation  as  “the  energy  lost  upon 
transmission  of  a  compressional  wave  from  all  causes  and  is  thus 
‘effective  attenuation’  in  most  cases  below  the  top  few  meters  of 
sediments.”  Since  3.5  kHz  will  not  penetrate  deeply  into  the  sand 
layer,  it  may  be  prudent  to  use  an  estimate  of  attenuation 
characteristic  of  the  upper  few  meters  of  the  in  situ,  saturated 
sediment.  The  attenuation  for  just  the  top  few  meters  of  the  sand 
would  be  higher  due  to  interstitial  sea  water.  However,  this  study 
will  use  an  attenuation  value  of  0.5  dB/m-kHz.  Hamilton’s  (1980) 
attenuation  (kp)  versus  depth  plot  is  given  in  Figure  3  to  illustrate 
the  depth  dependence  of  attenuation  for  sand  and  silt-clay 
sediments.  Notice  that  initially  the  silt-clay  attenuation  increases 
with  depth  and  then  decreases  while  the  sand  sediment  shows  a 
consistent  decrease  in  attenuation  with  depth. 

Shear  wave  speed  and  shear  attenuation  were  not  calculated 
for  this  study  since  shear  wave  speeds  are  very  low  relative  to 
compressional  wave  speeds  and  shear  attenuation  is  presumed  to  be 
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negligible.  Also,  the  transmission  loss  models  employed  do  not  use 
these  parameters  as  input  variables. 

The  complete  NPS  geoacoustic  model  for  the  sediment  layer  is 
given  in  Table  3  and  the  compressional  wave  speed  versus  depth  is 
plotted  in  Figure  10  which  incorporates  the  water  sound  speed 
profile  measured  for  Run  18. 

D.  LAMBERT’S  LAW  OF  BOTTOM  BACKSCATTERING 

Lambert’s  law  relationship  between  backscattering  strength 
and  grazing  angle  has  been  found  to  be  a  good  approximation  to 
measured  data  for  many  deep  water,  very  rough  bottoms  at  grazing 
angles  below  about  45°  (Urick,  1983).  Moreover,  this  relationship 
between  backscattering  strength  and  angle  appears  to  be  frequency 
independent.  Mackenzie’s  (1961)  optimum  value  for  “k”  (i.e.,  10  log 
p  =  k)  was  about  -29  dB  which  he  felt  adequately  described  the 

scattering  strength  of  clay,  mud,  silt  and  fine  sand  bottoms  at 
frequencies  between  0.5  kHz  and  100  kHz.  However,  the  predictable 
behavior  given  in  Lambert’s  law  applies  only  to  what  Urick  (1983) 
refers  to  as  a  Type  III  acoustic  bottom,  i.e.,  a  heavily  dissected 
bottom  with  underwater  ridges  for  which  Lambert’s  law  is  invariant 
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in  frequency.  However,  various  data  suggest  a  definite  frequency 
dependency  for  smoother  (mud,  clay  and  sand)  bottoms  (Urick,  1983). 
Sand  and  silt  bottoms  show  a  rise  in  scattering  strength  with 
frequency  at  the  rate  of  about  3  dB/octave  (or  a  variation  as  the 
first  power  of  the  frequency),  whereas  little  or  no  frequency 
dependency  appears  to  exist  for  rock,  sand  and  rock,  or  silt  and  shell 
bottoms,  Urick  attributes  this  difference  in  frequency  dependency 
as  possibly  due  to  the  difference  in  the  scale  of  bottom  roughness. 
Bottoms  with  roughnesses  large  compared  to  a  wavelength  have  a 
backscattering  strength  independent  of  frequency.  Bottoms  with  an 
appreciable  portion  of  their  roughness  spectrum  at  roughnesses  less 
than  an  acoustic  wavelength  will  have  a  scattering  strength  which 
increases  with  frequency.  Later  in  this  study,  measured 
reverberation  data  will  be  compared  to  modeled  reverberation  levels 
using  various  values  for  10  log  p,  to  determine  if  bottom  scattering 
strength  follows  Lambert’s  law  in  Area  Foxtrot. 
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Table  1.  The  current  geoacoustic  model  parameters  for  Area  Foxtrot. 


Compressional  Sound 
Speed  of  the  Sediment 

Density  of  the 
Sediment 

Attenuation  of  the 
Sediment 

1554  m/s 

2.0  g/cm3 

0.422  dB/m-kHz 
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Table  2.  Some  NPS  geoacoustic  model  parameters  for  the  sand 
bottom  in  Area  Foxtrot. 


Compressional 
Sound  Speed  of 
the  Sediment 

Density  of 
the  Sediment 

Mean  Grain  Size 
of  the  Sediment 
in  Phi  Units 

Attenuation 
of  the 

Sediment 

1817  m/s 

2.6  g/cm3 

1.33  phi  units 

0.5  dB/m-kHz 
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Table  3.  Compressional  sound  speed  versus  depth  in  the  NPS 
geoacoustic  model  -  includes  water  column  and  sediment  layer. 


Depth 

Speed 

Depth 

Speed 

(m) 

(m/s) 

(m) 

(m/s) 

0 

1519.00 

73 

1937.78 

6.25 

1518.00 

74 

1938.18 

12.5 

1516.00 

75 

1938.57 

18.75 

1512.00 

76 

1938.95 

25 

1505.00 

77 

1939.33 

31.25 

1506.00 

78 

1939.71 

37.5 

1500.00 

79 

1940.08 

43.75 

1496.00 

80 

1940.44 

50 

1493.00 

81 

1940.81 

56.25 

1490.00 

81 

1940.81 

62.5 

1488.00 

83 

1941.52 

64 

1484.00 

84 

1941.87 

64.0001 

1817.00 

85 

1942.21 

65 

1934.41 

86 

1942.55 

66 

1934.85 

87 

1942.89 

67 

1935.29 

88 

1943.22 

68 

1935.72 

89 

1943.55 

69 

1936.14 

90 

1943.88 

70 

1936.56 

91 

1944.20 

71 

1936.97 

92 

1944.52 

72 

1937.38 

93 

1944.83 

94 

1945.14 
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IV.  TRANSMISSION  LOSS  MODEL  SENSITIVITY  to  INPUT 
GEOACOUSTIC  PARAMETERS 


Now  that  a  detailed  Hamilton  geoacoustic  model  (the  NPS 
model)  has  been  developed  for  Area  Foxtrot,  a  TL  model  can  be  used 
to  accurately  calculate  TL  using  the  NPS  geoacoustic  parameters  as 
input.  The  TL  model  that  we  chose  to  evaluate  in  this  study  is  the 
GSM/MPE  model  provided  to  NPS  by  Dr.  Judith  Bishop  of  NUWC  with 
permission  from  Mr.  Ed  Chaika  of  AEAS.  GSM/MPE  is  a  relatively  new 
model  but  was  chosen  because  it  is  very  user-friendly  and  readily 
calculates  TL,  RL  and  other  active  sonar  model  parameters. 

A.  THE  GSM/MPE  PROPAGATION  MODEL 

The  GSM/MPE  (or  GSM)  model  allows  the  user  to  input 
variables  into  a  variety  of  sub-models  which  are  used  to  calculate 
TL  and  RL.  The  model  is  menu-driven  and  easy  to  use.  The  sub-model 
categories  are:  (1)  surface  reflection,  (2)  volume  attenuation,  (3) 
bottom  reflection,  (4)  bottom  phase  (5)  transmit  beam  pattern,  and 
(6)  receive  beam  pattern.  Each  of  the  six  sub-model  categories 
consists  of  one  or  more  sub-models  or  a  user-defined  function  which 
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is  selected  by  the  user.  User  inputs  include  minimum  and  maximum 
ranges,  range  increments,  minimum  and  maximum  frequency, 
frequency  increment  (if  a  frequency  range  is  used),  surface  and 
bottom  depths,  source  and  target  depths,  eigenray  tolerance, 
sediment  attenuation,  ocean/sediment  sound  speed  ratio, 
sediment/ocean  density  ratio,  wind  speed,  beam  patterns  and  ocean, 
sound  speed  profile.  Other  inputs  may  be  required  depending  on  the 
sub-models  selected. 

GSM  defaults  to  the  Rayleigh  reflection  coefficient  for  bottom 
loss  but  a  bottom  loss  versus  grazing  angle  curve  may  be  inserted  by 
the  user.  GSM  does  not  include  a  complete  geoacoustic  model  of  the 
sub-bottom  (sediment)  in  its  TL  calculations.  Only  geoacoustic 
parameters  at  the  surface  of  the  sediment  are  used  and  the  change  in 
these  parameters  with  depth  is  not  addressed.  This  may  be  a 
significant  weakness  in  using  GSM  in  shallow  water  since  it  will  be 
shown  later  in  this  study  that  the  sediment  layer  can  refract  a 
significant  amount  of  energy  back  into  the  water  column.  In  deep 
water  the  erroneous  concept  of  “negative  bottom  loss”  was 
introduced  because  the  impact  of  sub-bottom  energy  returned  to  the 
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water  column  was  not  addressed.  In  shallow  water  the  impact  of 
sub-bottom  returned  energy  should  be  greater  than  in  deep  water. 
Therefore,  this  returned  energy  must  be  modeled  accurately. 

Two  TL  calculations  were  performed,  one  for  the  current 
geoacoustic  model  (essentially  a  silt-clay  bottom)  (Table  1)  and  one 
for  the  NPS  geoacoustic  model  (sand  bottom)  (Table  2)  based  on 
geologic  data  acquired  near  Run  18.  Only  the  geoacoustic  parameters 
of  sediment  attenuation,  ocean/sediment  sound  speed  ratio  and 
sediment/ocean  density  ratio  in  the  bottom  reflection  (Rayleigh 
reflection)  sub-model  were  varied  between  the  two  TL  runs.  The 
variables  and  sub-models  used  for  each  TL  run  are  given  with  each 
TL  plot.  Note  that  a  source  depth  of  7.3  m  and  a  receiver  depth  of  64 
m  (same  as  bottom  depth)  was  used  because  we  are  primarily 
concerned  with  bottom  backscattering  in  this  study.  The  placement 
of  the  receiver  on  the  ocean  floor  will  give  TL  calculated  from  the 
source  to  the  bottom  where  bottom  backscatter  occurs  (Figure  11). 

Figures  12  and  13  are  the  GSM  TL  plots  for  the  NPS  and  the 
current  (NUWC)  geoacoustic  models,  respectively.  Dr.  Henry 
Weinberg  of  NUWC  also  ran  these  TL  cases  to  verify  the  results. 
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Close  comparison  of  Figures  12  and  13  shows  that,  as  expected, 
poorer  propagation  occurs  for  the  current  geoacoustic  model.  At  15 
km  the  NFS  model  shows  about  75  dB  loss  where  the  current 
geoacoustic  model  shows  about  81  dB  loss,  a  difference  of  about  6 
dB.  At  20  km  the  current  geoacoustic  model  shows  about  8  dB  more 
loss  than  the  NFS  model.  These  plots  are  presented  to  demonstrate 
the  high  sensitivity  of  TL  between  the  two  geoacoustic  models  and 
shows  that  the  silt-clay  sediment  yields  significantly  more  loss 
than  the  sand  sediment.  Recall,  Area  Foxtrot  sediments  are 
spatially  variable  and  no  one  geoacoustic  model  can  be  applied  to  the 
entire  area. 

The  above  calculations  were  performed  using  an  eigenray 
tolerance  of  0.10.  Eigenray  tolerance  reduces  (increases)  the 
running  time  of  the  multipath  eigenray  model  by  increasing 
(decreasing)  the  convergence  tolerance.  Convergence  is  assumed 
when  the  difference  in  random  addition  between  consecutive  ray 
cycles  is  less  than  the  eigenray  tolerance  indicated  (Weinberg, 

1985).  The  TL  calculations  were  performed  again  but  with  an 
eigenray  tolerance  of  0.0010  to  try  and  increase  the  accuracy  of  the 
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TL  by  including  more  eigenrays  in  the  calculations.  Figures  14  and 
15  are  the  GSM  TL  plots  with  the  lower  eigenray  tolerance  for  the 
NFS  and  current  geoacoustic  models,  respectively.  Comparison  of 
Figures  12  and  14  (NFS  model)  and  of  Figures  13  and  15  (current 
model)  shows  very  little  difference  in  TL  between  the  two  eigenray 
tolerance  settings.  However,  since  the  calculation  time  is  only  a 
few  seconds  greater  for  the  smaller  eigenray  tolerance,  it  was  used 
for  all  subsequent  TL  calculations  in  this  thesis. 

The  TL  was  then  calculated  for  both  the  NFS  and  current 
geoacoustic  models  using  the  Finite  Element  Farabolic  Equation 
(FEFE)  model  to  illustrate  the  sensitivity  of  TL  to  bottom 
geoacoustic  properties.  The  results  are  discussed  in  the  next 
section. 

B.  THE  FEPE  PROPAGATION  MODEL 

The  FEPE  model  has  been  used  at  NFS  and  elsewhere  for  several 
years  and  has  been  found  to  be  very  accurate  by  a  large  number  of 
scientists  (Chin-Bing  et  al.,  1993).  It  is  less  user-friendly  than  the 
GSM  model  and  takes  much  longer  to  run  but  has  the  advantage  of 
including  a  full  Hamilton-type  geoacoustic  model  of  the  sub-bottom 
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as  input.  That  is,  where  the  GSM  model  uses  only  the  compressional 
wave  speed  and  density  ratios  across  the  water/sediment  interface, 
and  one  value  of  sediment  attenuation,  FEPE  allows  the  sediment 
compressional  wave  speed,  density,  and  attenuation  to  vary  with 
depth  and  range.  As  such,  FEPE  accounts  for  modes  propagated  into 
and  out  of  the  sediment  layer.  Although  FEPE  is  normally  used  with 
lower  frequencies,  it  can  be  used  for  our  frequency  of  interest,  3.5 
kHz  (Campbell,  1995).  Care  must  be  taken  however  when  selecting 
the  resolution  of  the  calculations,  in  this  study,  a  range  increment 
of  one  quarter  acoustic  wavelength  (A,/4  )  and  a  depth  resolution  of 

one  eighth  acoustic  wavelength  (?i/8)  was  used  (Campbell,  1995). 

FEPE  can  also  be  used  in  range-dependent  seafloor  environments. 

The  bottom  in  Area  Foxtrot  is  range-dependent  (Figures  6,  7,  8,  and 
9).  However,  since  the  purpose  of  these  TL  estimates  is  to  show  the 
impact  of  bottom  geoacoustic  properties  on  TL  (and  ultimately  RL),  a 
range-independent  bottom  environment  is  assumed. 

The  NPS  model  input  to  FEPE  (Table  3)  incorporated  a  depth- 
dependent  compressional  wave  speed,  density  and  attenuation  in  the 
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30  m  sediment  layer  according  to  Hamilton’s  equations  (Hamilton, 
1976,  1979).  Below  the  sediment  layer  is  an  attenuating  half-space. 
The  current  geoacoustic  model  input  to  FEPE  was  the  same  as  that 
input  to  GSM.  For  this  situation,  no  depth  dependence  of  geoacoustic 
parameters  were  permitted  since  they  were  not  part  of  the  current 
model  (and  are  not  considered  in  GSM). 

Figure  16  shows  the  TL  results  from  FEPE  for  the  NPS  and 
current  geoacoustic  models  for  a  source  at  7.3  m  and  receiver  at  64 
m.  Notice  the  significant  difference  in  TL  as  range  increases.  At  5 
km  the  current  model  shows  about  10  dB  more  loss  than  the  NPS 
model,  at  10  km  and  15  km  the  difference  increases  to  about  18  dB 
and  23  dB,  respectively.  For  comparison,  TL  was  calculated  using 
FEPE  for  both  geoacoustic  models  but  with  the  source  and  receiver 
both  at  7.3  m  (Figure  17).  The  TL  for  each  geoacoustic  model  is 
slightly  less  than  that  when  the  receiver  is  at  the  bottom  (64  m). 
However,  the  relative  difference  between  TL  for  the  two  geoacoustic 
models  is  the  same  whether  the  receiver  is  at  the  bottom  or  at  the 
source  depth.  Figure  18  is  a  combination  of  Figures  14  and  15  and 
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shows  TL  for  the  NPS  and  current  geoacoustic  models  for  each 
source/receiver  geometry. 

This  dramatic  difference  in  TL  with  range  when  modeling  a 
“hard”,  coarse  sand  bottom  with  a  sediment  compressional  wave 
speed  of  1817  m/s  (NPS  model)  and  a  “soft”  bottom  with  a  sediment 
compressional  wave  speed  of  1554  m/s  (current  model)  shows  the 
importance  of  knowing  the  correct  geology  and  geophysics  of  a  study 
area  before  constructing  geoacoustic  models.  Figures  6  through  9 
show  the  spatial  variability  of  surface  sediment  type  in  Area 
Foxtrot.  The  sediment  type  can  also  vary  with  depth  although  these 
data  were  not  readily  available  for  this  study. 

C.  COMPARISON  OF  GSM  AND  FEPE  TL  CALCULATIONS  FOR  THE 
SAME  GEOACOUSTIC  MODEL 

TL  curves  calculated  using  the  current  geoacoustic  model  as 
input  to  the  GSM  and  FEPE  models  were  plotted  together  for 
comparison  (Figures  19  and  20)  in  order  to  quantify  and  illustrate 
the  difference,  if  any,  in  their  TL  estimates.  The  GSM  TL  energy 
levels  were  summed  both  coherently  (phased  summation)  and 
incoherently  (random  summation)  while  the  FEPE  levels  were 
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summed  coherently.  These  figures  show  that  the  GSM  and  FEPE  TL 
curves  are  in  agreement  out  to  about  2.5  km  in  range.  Beyond  that, 
the  curves  diverge  and  show  that  the  GSM  TL  curve  underestimates 
the  loss  relative  to  FEPE  which  incudes  the  full  geoacoustic  model. 
At  a  range  of  10  km  the  GSM  TL  is  about  7  dB  less  than  the  FEPE  TL 
and  at  15  km  the  difference  has  increased  to  about  13  dB  less  loss. 

The  difference  in  the  TL  calculations  and  divergence  of  the  TL 
curves  is  likely  due  to  how  the  models  handle  the  physics  of 
reflection  at  the  bottom  interface  and  refraction  into  and  out  of  the 
sub-bottom  sediments  as  well  as  the  modeled  geology.  The  soft 
bottom  current  model  has  a  critical  angle  at  the  sediment  interface 
of  17°  which  means  that  more  energy  from  an  omnidirectional 
source  will  penetrate  into  the  bottom  than  for  a  harder  bottom 
which  will  have  a  larger  critical  angle.  The  two  TL  curves  agree  at 
short  ranges  (less  than  2.5  km),  for  high  grazing  angles  and  for 
relatively  few  bottom  bounces.  But  as  range  and  the  number  of 
bottom  bounces  increase,  GSM  does  not  account  for  enough  energy 
loss  into  the  sediment  layer.  Because  GSM  underestimates  the  loss 
per  bottom  interaction  and  overestimates  the  reflection  coefficient 
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at  the  sediment  interface,  the  GSM  and  FEPE  curves  diverge.  The 
term  “bottom  interaction”  is  applied  instead  of  “bounce”  because 
with  a  soft,  mud  bottom  energy  will  penetrate  the  sediment  and  be 
scattered  and  attenuated  within  the  sediment  during  refraction  back 
up  toward  the  water  column.  Hence,  not  all  the  energy  incident  on 
the  bottom  is  truly  “bounced”  from  the  interface.  Therefore,  as  the 
nurnber  of  bottom  interactions  increases,  the  separation  between 
the  TL  curves  increases.  Because  the  GSM  model  only  considers 
bottom  loss  from  the  water  /  sediment  interface  as  a  function  of 
grazing  angle  and  does  not  include  the  entire  geoacoustic  model  (i.e., 
within  the  sediment  layer),  it  fails  to  account  for  the  additional  TL 
due  to  attenuation  and  refraction  within  the  sediment.  FEPE  does 
account  for  the  entire  geoacoustic  model  of  the  sub-bottom  and 
attenuation  of  energy  refracted  through  the  sediment  and  back  into 
the  water  column. 

TL  estimates  calculated  for  the  NPS  geoacoustic  model  (hard 
sand  bottom)  using  GSM  (phased  and  random  TL)  and  FEPE  (phased  TL 
only)  are  compared  in  Figures  21  and  22.  Here,  the  TL  curves  show  a 
nearly  constant  difference  of  about  5  dB  out  to  15  km  range.  The 
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GSM  TL  estimates  show  more  loss  than  the  FEPE  TL  estimates  for 


the  NPS  geoacoustic  model  (hard  bottom).  Also,  the  two  TL  curves 
do  not  diverge  with  range  due  to  the  harder  bottom  modeled  where 
much  less  energy  is  refracted  through  the  sediment  layer  and  back 
into  the  water  column.  Therefore,  attenuation  within  the  sediment 
is  not  as  significant  a  factor  for  the  hard  bottom  and  bottom 
interactions  can  be  treated  as  “bottom  bounces”. 

Since  the  sub-models  within  GSM  can  not  accurately  estimate 
TL  in  shallow  water,  the  TL  curves  from  GSM  will  be  corrected  by 
adding  in  TL  corrections  determined  from  the  FEPE  estimates.  This 
will  allow  us  to  still  use  GSM  to  estimate  reverberation  levels  from 
the  corrected  TL  estimates. 

NUWC  has  addressed  this  problem  by  calculating  bottom  loss 
using  the  Seismo-Acoustic  Fast  Field  Algorithm  for  Range 
Independent  Environments  (SAFARI)  model  and  importing  these 
curves  into  GSM  to  estimate  TL  (Jensen,  1995).  This  procedure 
shows  some  merit  as  TL  estimates  from  GSM  using  the  SAFARI 
bottom  loss  curves  do  compare  well  to  some  measured  TL  curves  for 
a  small  part  of  Run  18.  The  FEPE  TL  curves  also  compare  well  to  the 
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measured  TL.  It  must  be  noted  that  the  TL  measurements  were 
collected  between  the  surface  ship  and  the  target  submarine  (Figure 
8).  Energy  propagation  through  this  path  is  over  both  sand  and  silt- 
clay  bottoms  (Figure  7).  Therefore,  more  than  one  geoacoustic  model 
is  required  to  properly  estimate  TL. 

NFS  has  chosen  to  use  FEPE  (instead  of  SAFARI)  to  augment 
GSM  because  it  accounts  for  energy  that  is  refracted  within  the 
sediment  and  returned  to  the  water  column.  FEPE  is  also  likely  to 
become  the  Navy  standard  Oceanographic  and  Atmospheric  Master 
Library  (OAML)  approved  shallow  water  TL  model.  However,  the  main 
point  stressed  by  Figures  16  through  20  is  that  GSM  requires  an 
additional  model  (whether  it  be  SAFARI  or  FEPE)  to  correctly 
account  for  the  sub-bottom  interacting  energy,  even  at  3.5  kHz.  FEPE 
correctly  addresses  the  propagation  of  the  sub-bottom  refracted 
energy  that  is  returned  to  the  water  column.  This  avoids  the  need  to 
incorporate  the  “negative  bottom  loss”  energy  concept  introduced  in 
the  mid-1980’s  to  overcome  the  fact  that  energy  refracted  back  into 
the  water  column  by  sediments  was  neglected  in  deep  water  bottom 
loss  models. 
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V.  BOTTOM  BACKSCATTER 


A.  THE  SAIC  METHOD  OF  ANALYSIS  TO  DETERMINE  BOTTOM 
SCATTERING  STRENGTH  (AND  THE  DOMINANT  BOTTOM 
SCATTERING  MECHANISM) 

Estimation  of  bottom  scattering  strength  is  normally 
conducted  by  direct  measurement  methods  (Urick,  1983;  McCammon', 
1993).  Direct  measurements  of  bottom  backscattering  strength  can 
yield  misleading  results  if  the  acoustic  /  geologic  environment  of 
the  study  area  is  not  well  known  (this  point  is  addressed  in  detail 
later  In  this  study).  SAIC  has  developed  a  method  by  which  bottom 
scattering  strength  may  be  deduced  through  comparison  of  measured 
and  modeled  RL.  This  method  has  the  advantage  of  including  full 
Hamilton-type  geoacoustic  models  in  the  analysis. 

The  following  reverberation  data  analysis  technique  was 
conducted  by  SAIC  (Hanna  et  al.,  1994)  as  part  of  the  Harsh 
Environments  Program  (HEP)  of  the  Office  of  Naval  Research.  Their 
methods  and  results  are  described  below. 
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B.  DETERMINATION  OF  SEAFLOOR  SCATTERING  STRENGTH  AND 
BACKSCATTER 

The  first  point  to  establish  in  the  analysis  of  the  data  is  to 
determine  which  bottom  scattering  mechanism  controls  bottom 
reverberation  in  Area  Foxtrot,  realizing  that  more  than  one 
mechanism  may  be  present.  Scattering  may  be  either  at  the  water- 
bottom  interface  or  within  the  bottom  sediments.  Further,  if  the 
dominant  scattering  mechanism  is  found  to  be  within  the  bottom 
sediments,  the  bulk  of  the  scatter  could  then  be  due  to  either 
inhomogeneities  within  the  sediment  or  from  layering  (macro  or 
micro)  within  the  sediment  column.  For  clarity,  scattering  from  the 
water-bottom  interface  will  be  referred  to  as  bottom  interface 
scatter.  Scattering  from  within  the  bottom  sediments  will  be 
referred  to  as  bottom  volume  scatter.  The  bottom  scattering 
strength  equation  can  then  be  written  as 


Sb  —  A(Sbi)  +  B  (Sbv) 


(6) 


where  Sb  is  the  total  bottom  scattering  strength, 
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A  and  B  are  relative  percentages, 

Sbi  is  the  scattering  strength  at  the  ocean-sediment 
interface,  and 

Sbv  is  the  scattering  strength  within  the  sediment  volume. 

Due  to  the  nature  of  this  study,  we  assume  the  reverberation- 
limited  case  of  the  active  sonar  equation  to  be: 


SE  =  (SL  -  TLst  +  TS  -  TLtr)  -  RL  -  DT  (7) 

RL  =  (SL  -  TLsb  +  Sb  +  10  log(A)  -  TLbr)  (8) 

where  SE  is  the  signal  excess, 

SL  is  the  source  level, 

TLst  is  the  transmission  loss  along  source  to  target  path, 
TS  is  the  target  strength, 

TLtr  is  the  transmission  loss  along  target  to  receiver  path, 
RL  is  the  reverberation  level. 


65 


TLsb  is  the  transmission  loss  from  source  to  bottom  path, 
Sb  is  the  scattering  strength, 

A  is  the  reverberating  bottom  area 

TLbr  is  the  transmission  loss  from  the  bottom  scattering 
area  along  the  path  to  receiver,  and 
DT  is  the  detection  threshold. 

Equation  7  is  often  written  in  a  specialized  form  for  shallow 
water  cases.  In  this  case,  it  is  assumed  that  the  reverberation 
originates  from  the  bottom  interface  and  that  the  propagation  path 
to  that  bottom  area  is  the  same  as  the  path  to  the  target  at  equal 
range.  Using  this  assumption,  all  the  propagation  terms  in  Equations 
7  and  8  cancel  leaving  signal  excess  a  function  of  target  strength, 
detection  threshold,  scattering  strength  and  the  area  term.  Equation 
7  then  becomes: 


SE  =  TS  -  DT  -  Sb  -  10  log(A) 


(9) 


and  the  reverberation  equation  is  then 


RL  =  SL  -  2TL  +  Sb  - 10  log(A)  (1 0) 


where  2TL  is  the  2-way  transmission  loss. 


Here,  we  assumed  that  TLst  =  TLsb,  and  TLtr  «  TLbr-  Implicit  in  these 
approximations  is  that  the  scattering  mechanism  is  of  the  bottom 
interface  type.  Earlier  studies  used  these  same  assumptions  in 
attempting  to  estimate  the  bottom  scattering  strength  (Urick, 

1983).  Note  that  since  the  direction  of  scatter  is  back  to  the  source 
for  monostatic  sonars,  scattering  strength  can  be  more  specifically 
referred  to  as  backscattering  strength. 

The  SAIC  method  starts  by  fitting  the  measured  propagation 
loss  from  the  source  to  the  receiver  to  a  curve.  The  reverberation 
level  s  then  calculated  using  existing  models  under  the  assumption 
of  zero  scattering  strength  (Sb  =  0,  or  SS=0  in  Hanna  et  al.,  1994 
terminology).  The  measured  propagation  curve  (data)  is  then 
subtracted  from  the  modeled  reverberation  curve  (Sb  =  0).  The 
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difference  between  the  two  curves  (Modeled  -  Data)  is  assumed  to  be 
the  backscattering  strength  if  all  other  reflection  and  scattering 
mechanisms  are  properly  accounted  for  in  the  model.  As  seen  in 
Figure  23,  the  resulting  backscattering  strength  curve  is  range 
dependent  because  the  modeled  and  measured  decay  rates  are 
obviously  different.  However,  this  implies  that  the  scattering 
strength  is  not  only  a  function  of  backscattering  at  the  bottom 
interface,  but  also  may  include  sub-bottom  volume  scattering 
geoacoustics  (Mourad  and  Jackson,  1993;  Stephen,  1993;  Ellis  et  al., 
1993;  Jackson,  1992).  Scattering  strength  should  not  have  a 
position  dependence  that  varies  monotonically  with  distance  to  the 
receiver  assuming  a  homogeneous  bottom.  The  reverberation  data  in 
Figure  23  are  for  the  octave  band  centered  at  100  Hz  and  are  for  the 
Nile  Fan,  but  similar  results  have  been  found  to  hold  for  the  Korean 
shallow  water  areas  (Hanna  et  al.,  1994). 

Figure  24  shows  the  same  analysis  as  the  previous  figure  but 
for  800  Hz.  This  implies  that  the  scattering  strength  anomaly  is 
seen  across  the  frequency  band  from  100  Hz  to  800  Hz. 
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For  sediment  volume  scattering  to  be  important  in  the 
scattering  strength  term,  energy  must  penetrate  into  the  ocean 
bottom.  In  ray  theory  terminology,  this  implies  that  acoustic  energy 
must  penetrate  the  bottom  interface  with  an  angle  greater  than  the 
critical  angle.  In  shallow  water  regions,  the  sediment  is  often 
composed  of  at  least  some  fractional  portion  of  sand  which  causes 
the  sound  speed  ratio  at  the  ocean-sediment  interface 
(Cwater/Csediment  =  n)  to  be  less  than  unity  and  the  ratio  of  densities 
(Psediment  /  Pwater  =  Hi)  to  be  greater  than  n.  This  case  is  most 
commonly  found  in  naturally  deposited  bottom  sediments  and  results 
in  complete  reflection  from  the  interface  for  grazing  angles  less 
than  the  critical  angle  (Urick,  1983  p.  137).  In  the  case  of  a  bottom 
comprised  of  a  soft  mud  or  saturated  ooze,  the  Cwater/Csediment  ratio 
may  be  greater  than  unity  and  little  reflection  takes  place;  most  of 
the  energy  is  transmitted  into  the  sediment  until  an  angle  of 
intromission  is  reached  when  all  of  the  energy  is  transmitted  into 
the  bottom. 

Thus,  while  bottom  interface  scattering  can  effect  all  rays 
impinging  on  the  bottom,  bottom  volume  scatter  can  only  effect 
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those  rays  which  have  bottom  grazing  angles  greater  than  the 
critical  angle.  Figure  25  shows  that  for  bottom  interacting  rays, 
energy  loss  associated  with  angles  below  the  critical  angle  will 
decay  with  range  less  rapidly  than  that  for  rays  greater  than  the 
critical  angle.  This  difference  results  because  acoustic  energy  that 
penetrates  into  the  sediment  layer  will  be  attenuated  at  a  faster 
rate  due  to  the  higher  attenuation  coefficient  of  the  sediment 
compared  to  the  water. 

Consequently,  we  will  take  measured  reverberation  levels  and 
subtract  out  modeled  reverberation  levels  for  which  the  scattering 
strength  (Sb)  is  assumed  to  be  zero.  What  remains  should  be  the 
actual  bottom  scattering  strength  (Sbi  +  Sbv)-  By  inspection  of  the 
energy  loss  decay  as  a  function  of  range  in  the  bottom  scatter 
curves,  we  will  be  able  to  determine  whether  bottom  interface  or 
bottom  volume  scatter  is  the  dominant  scattering  mechanism. 

Having  determined  the  dominant  scattering  mechanism  at  the 
bottom,  we  can  better  evaluate  the  bottom  scattering  kernel  and 
physics  involved  in  backscatter  from  the  bottom. 
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C.  USING  THE  SAIC  ANALYSIS  METHOD  TO  DETERMINE  BOTTOM 
SCATTERING  STRENGTH  (Sb)  FROM  AREA  FOXTROT  DATA 

The  GSM  model  was  run  to  calculate  the  reverberation  level 
using  the  NPS  geoacoustic  model  with  the  bottom  scattering 
strength  set  to  zero  as  described  above.  The  model  results  were 
then  subtracted  from  the  measured  RL  data.  Theoretically,  the 
remaining  curve  should  describe  the  bottom  scattering  strength  as  a 
function  of  time.  Figures  26,  27  and  28  compare  the  modeled  RL 
(Sb  =  0)  with  the  measured  RL  data  (upper  panels)  and  the  difference 
between  the  two  curves  for  beams  1,  12  and  24,  respectively,  for 
ping  2  of  Run  18  (lower  panels).  Also  shown  is  a  3rd-order 
polynomial  fit  to  the  reverberation  level  difference.  This  curve  can 
be  used  as  a  measure  of  the  bottom  scattering  strength.  Note  the 
difference  between  the  modeled  and  measured  reverberation  levels 
for  each  beam.  The  plots  for  beams  1  and  24  show  about  a  65  dB 
difference  at  20  seconds  whereas  beam  12  shows  about  a  70  dB 
difference.  Also,  the  bottom  scattering  strength  curve  for  beam  12 
has  a  greater  slope  than  the  curves  for  beams  1  and  24.  This 
indicates  that  the  bottom  scattering  strength  for  beam  12  more 
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closely  approximates  that  given  by  Lambert’s  law  than  the 
scattering  strength  for  the  other  two  beams.  Recall  that  beams  1 
and  24  propagate  mostly  over  a  silt-clay  bottom  whereas  beam  12 
was  mostly  over  a  coarse  sand  bottom  (Figures  6  through  9).  This 
difference  in  sediment  properties  causes  an  azimuthal  dependence 
of  bottom  scattering  strength  which  in  turn  causes  an  azimuthal 
dependence  in  the  reverberation  level. 

The  RL  was  next  calculated  using  the  NFS  sand  bottom  model 
as  input  with  the  bottom  scattering  strength  described  by  a 
Lambert’s  law  formulation: 

Sb  =  -1 9.5  +  10  log  sin^  0 .  (11) 


These  RL  values  were  compared  with  measured  RL  data  for  the  odd- 
numbered  beams  for  pings  2,  30  and  60  (Figures  29  (a-d),  30  (a-d) 
and  31  (a-d),  respectively),  to  see  if  the  azimuthal  dependence  of 
reverberation  was  a  consistent  feature.  As  seen  in  these  figures, 
azimuthal  dependence  of  reverberation  is  evident  from  ping  2 
through  to  ping  60.  Note  that  the  GSM  modeled  RL  is  not  very 
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accurate  due  to  the  inherent  limitation  of  the  model  in  estimating 
TL  as  discussed  earlier. 

To  further  examine  the  dependence  of  reverberation  level  on 
azimuth  angle,  the  measured  reverberation  level  was  plotted  as  a 
reverberation  surface  (ignoring  beam  angle  in  the  display)  for  all 
beams  for  pings  2,  20,  40,  60  and  88  (Figures  32,  33,  34,  35  and  36, 
respectively).  This  essentially  shows  how  RL  changed  as  a  function 
of  distance,  i.e.,  as  the  acoustic  source  moved  eastward  during  Run 
18.  Figure  37  is  a  three-dimensional  surface  plot  of  the 
reverberation  field  for  ping  88,  approximately  5  km  distance  from 
ping  2.  These  reverberation  surfaces  show  the  azimuthal  dependence 
of  reverberation  very  well.  It  is  easily  seen  in  Figures  32  through 
37  that  as  time  increases,  reverberation  levels  for  the  beams 
steered  ahead  towards  the  harder,  sand  bottom  (beams  6-19) 
maintain  higher  levels  than  the  peripheral  beams  (beams  1-5  and 
20-24)  steered  toward  the  softer,  mud  bottom.  Beams  to  the  right 
of  center  (beams  1-5)  show  the  greatest  reduction  in  reverberation 
with  time  and  also  that  the  degree  of  change  across  the  beam 
spectrum  changes  from  ping  to  ping.  This  implies  an  inhomogeneous 
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bottom  in  the  area  of  Run  18  which  is  consistent  with  the  NPS 
geoacoustic  contours  in  Figure  7. 

The  azimuthal  dependence  of  reverberation  can  be  explained  by 
the  local  geology.  Referring  back  to  the  surface  sediment 
compressional  wave  speed  plots  (Figures  7  and  8)  it  is  seen  that  for 
Run  18  the  more  central  beams  of  the  sonar  interact  with  a  higher 
speed  bottom  for  a  much  longer  time  than  the  peripheral  beams.  The 
peripheral  beams  interact  with  a  slower,  silt-to-clay  bottom  which 
lies  to  the  left  and  right  of  the  high  speed  sediment.  This  high  speed 
bottom  forms  essentially  a  strip  which  lies  somewhat  along  the 
track  of  the  surface  ship.  The  slowest  speed  bottom  is  found  to  the 
right  of  the  surface  ship  track  where  the  lower  number  beams  (1-5) 
interact  with  the  bottom.  This  increased  bottom  penetration  for  the 
peripheral  beams  explains  the  reduction  in  reverberation  level  with 
time  and  also  explains  why  the  reverberation  trend  for  the 
peripheral  beams  does  not  follow  Lambert’s  law.  Energy  which 
impinges  the  soft  bottom  will  enter  the  sediment  and  be  attenuated 
within  the  sediment.  Bottom  scatter  in  this  case  will  not  obey 
Lambert’s  law  and  is  assumed  to  be  a  function  of  sin  0  (Mourad  and 
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Jackson,  1993;  Jackson  and  Briggs,  1992)  vice  sin2  e.  The  value  of 


10  log  |j,  may  also  vary  and  is  in  part  a  function  of  the  interface 
roughness. 

To  quantify  the  bottom  scattering  strength,  reverberation 
levels  from  GSM  need  to  be  adjusted  to  account  for  the  inherent 
errors  in  the  GSM  TL  estimates  (when  bottom  loss  is  calculated  from 
Rayleigh  reflection  coefficients).  This  is  most  easily  done  by 
measuring  the  difference  in  TL  between  the  GSM  and  FEPE  estimates, 
then  multiplying  this  difference  by  two  (to  account  for  two-way 
travel  time)  and  adding  this  quantity  to  the  GSM  reverberation 
levels.  This  method  of  correcting  the  GSM  reverberation  levels  is 
consistent  with  Equation  10.  However,  due  to  the  lateral 
inhomogeneity  of  the  surface  sediment  (and  thus  compressional 
speed),  a  second  Hamilton-type  geoacoustic  model  was  developed  in 
addition  to  the  NPS  sand  bottom.  This  model  is  for  a  silt-clay 
bottom  with  a  surface  sediment  compressional  wave  speed  of  1580 
m/s.  The  FEPE  model  was  also  run  to  estimate  the  TL  for  the  NPS 
silt-clay  bottom  which  included  the  full  Hamilton  geoacoustic 
model. 
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GSM  reverberation  levels  were  corrected  (for  TL)  for  beams  1, 
12,  and  24  of  pings  2,  40,  60,  and  88.  These  beams  were  chosen 
because  they  were  likely  to  be  steered  in  the  least  range-dependent 
azimuths  (Figure  7).  Different  bottom  scattering  strength  kernels 
were  incorporated  into  the  NFS  sand  (beam  12)  and  NFS  silt-clay 
(beams  1  and  24)  geoacoustic  models.  TL  was  corrected  as 
described  earlier.  Beams  1  and  24  were  modeled  for  propagation 
over  the  NFS  silt-clay  bottom  while  beam  12  was  propagated  over 
the  NFS  sand  bottom.  Figures  38  (a-e),  39  (a-e),  and  40  (a-e)  show 
measured  reverberation  levels  plotted  with  uncorrected  GSM  RL  in 
the  upper  panel,  and  measured  reverberation  levels  plotted  with  GSM 
RL  with  corrected  TL  errors,  for  three  different  bottom  scattering 
strength  kernels  (bottom  panel).  On  initial  inspection  of  the  ping  2 
plots,  for  beam  1  (Figure  38a)  (NFS  silt-clay  bottom)  the  “best  fit” 
bottom  scattering  strength  kernel  appears  to  be  Sb  =  -35  +10  log  sin 

(0).  For  ping  2,  beam  12  (Figure  39a)  (NFS  sand  bottom)  it  appears 

that  none  of  the  proposed  bottom  scattering  strength  kernels  fits 
the  measured  data  throughout  the  entire  20  second  reverberation 
time  series.  For  ping  2,  beam  24  (Figure  40a)  (NFS  silt-clay 
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bottom)  the  “best  fit”  bottom  scattering  strength  kernel  appears  to 
be  Sb  =  -33  +10  log  sin  (0).  Similar  observations  are  made  for  pings 

60  and  88.  However,  close  examination  of  Figures  38a,  39a,  and  40a 
in  conjunction  with  the  best  geologic  and  geoacoustic  data  available 
(Figures  5  and  6)  show  some  very  interesting  results. 

For  Instance,  beams  1  and  24  (Figures  38a  through  38e  and  40a 
through  40e)  show  that  modeled  RL  for  a  silt-clay  bottom  with  a 
bottom  scattering  strength  kernel  employing  10  log  sin  (0)  generally 

fits  the  measured  RL  very  well  throughout  the  entire  time  series. 

This  kernel  is  consistent  with  a  bottom  volume  scattering 
mechanism  proposed  by  Jackson  and  Briggs  (1992)  and  Mourad  and 
Jackson  (1993).  But  Figures  6  and  7  show  the  bottom  in  the 
immediate  vicinity  of  the  surface  ship  to  be  the  NFS  sand  type 
bottom.  The  bottom  type  does  change  with  increased  range  becoming 
more  like  the  NFS  silt-clay  bottom  along  the  azimuths  of  beams  1 
and  24.  The  bottom  scattering  kernel  for  a  silt-clay  bottom  is  not 
expected  to  fit  the  entire  measured  RL  time  series  accurately  when 
the  initial  propagation  of  ping  2  is  over  hard  sand. 
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To  further  investigate  this  range-dependence  in  Sb,  the  GSM  RL, 
using  the  NFS  sand  bottom  as  input,  was  compared  with  the  GSM  RL 
for  the  silt-clay  bottom  and  with  the  measured  RL  for  beam  1  of 
ping  2  (Figure  41).  Figure  39  shows  that  for  beam  1  of  ping  2,  GSM 
RL  (with  corrected  TL)  for  the  NFS  sand  bottom  using  a  bottom 
scattering  strength  kernel  of  -33  +10  log  sin2  (0)  fits  the  measured 

RL  as  well  as  the  -35  +10  log  sin  (9)  kernel  out  to  about  7.5  seconds. 
From  7.5  seconds  to  20  seconds  the  -35  +10  log  sin  (0)  kernel  fits 
the  measured  RL  whereas  the  -33  +  10  log  sin2  (0)  kernel  over¬ 
estimates  RL.  For  beam  1  of  ping  2,  it  appears  that  at  least  two 
different  geoacoustic  models  (sand  and  silt-clay)  are  required  to 
describe  the  bottom  over  which  the  beam  propagates.  And  each 
geoacoustic  model  requires  a  different  bottom  scattering  strength 
kernel  to  estimate  RL.  For  the  sand  bottom  in  the  vicinity  of  the 
ping  origin  (out  to  7.5  sec)  -33  +10  log  sin2  (0)  appears  to  be  the 
appropriate  bottom  scattering  strength  kernel.  However,  for  the 
silt-clay  bottom  (beyond  7.5  sec)  -35  +10  log  sin  (0)  better 
describes  bottom  scattering  strength.  Using  these  two  geoacoustic 


78 


models  is  consistent  with  the  geology  of  Area  Foxtrot  (Figures  6  and 

7). 

There  are  three  variables  which  account  for  the  shape  of  the  RL 
curve  (RL  as  a  function  of  time  or  range)  in  Figure  41:  1)  the  TL  as  a 
function  of  renge,  2)  Mackenzie’s  coefficient  (10  log  jx),  and  3) 
whether  the  angular  dependence  is  sin  (0)  or  sin2  (e).  Factors  1  and 

3  above  may  be  linked  to  a  specific  sediment  type  (Jackson  and 
Briggs,  1992;  Mourad  and  Jackson,  1993).  A  silt-clay  bottom  should 
exhibit  a  TL  trend  that  increases  more  rapidly  with  range  (poorer 
propagation)  than  a  hard  sand  sediment  (Figures  16  and  17).  The 
silt-clay  sediment  is  also  a  possible  source  of  volume  reverberation 
(Jackson  and  Briggs,  1992;  Mourad  and  Jackson,  1993)  in  which  the 
bottom  scattering  strength  kernel  would  not  obey  Lambert’s  law 
(sin2  (0))  but  would  be  a  function  of  sin  (0).  This  is  due  to  the 

shallow  grazing  angles  and  “diffuse”  acoustic  impedance  across  the 
ocean  /  silt-clay  interface  which  allows  increased  energy 
penetration  into  the  sediment.  Note  that  an  angular  dependence  of 
sin  (0)  should  yield  a  higher  bottom  scattering  strength  than  for 
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sin2  (0),  but  the  sediment  attenuation  effectively  offsets  this 
increase  and  may  reduce  bottom  scattering  and  RL  to  levels  below 
those  found  when  the  angular  dependence  is  sin2  (0).  In  the  case  of 

Figure  41,  the  difference  between  the  10  log  |x  values  for  the  sand 
and  silt-clay  bottoms  is  2  dB  and  is  negligible.  This  suggests  that  it 
is  the  TL  trend  and  angular  dependence  that  dominates  the  RL  trend. 

For  comparison,  the  sand  bottom  was  then  assigned  a  bottom 
scattering  strength  kernel  of  -35  +10  log  sin  (0)  and  the  silt-clay 
bottom  was  assigned  a  bottom  scattering  strength  kernel  of  -33  +10 
log  sin2  (0),  i.e.,  the  bottom  scattering  strength  kernels  were 

switched  for  the  two  models.  This  separates  the  TL  trend  and  the 
angular  dependence  of  RL.  As  seen  in  Figure  42  neither  of  the  bottom 
type/bottom  scattering  strength  kernel  combinations  correctly 
model  the  rheasured  RL.  Remembering  that  there  are  three  variables 
presented  here  which  affect  the  RL  trend  (with  time  or  range), 

10  log  (j.  was  varied  in  the  bottom  scattering  strength  kernel  for  the 
silt-clay  bottom  and  employed  sin2  (0)  (Lambert’s  law)  in  the  kernel 
to  see  how  well  it  could  fit  GSM  RL  to  the  measured  RL. 
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Surprisingly,  GSM  RLs  for  the  silt-clay  bottom  using  two  very 
different  bottom  scattering  strength  kernels  (-35  +10  log  sin  (0)  and 

-17  +10  log  sin2  (e))  are  nearly  perfectly  coincident  (Figure  43). 

Figure  44  is  an  expanded  view  of  Figure  43  to  enhance  the  minor 
differences  in  the  RL  curves.  Figure  45  is  the  same  as  Figure  43 
except  the  sand  bottom  is  used  in  modeling  RL  with  the  bottom 
scattering  strength  kernel  of  -33  +10  log  sin2  (0)  (GSM  RL)  curve 
also  included. 

To  explain  the  similarity  of  the  RL  curves  in  Figure  43,  the  two 
bottom  scattering  strength  kernels  were  plotted  as  a  function  of 
grazing  angle  (Figures  46  and  47).  Notice  that  at  very  low  grazing 
angles  (less  than  about  3°)  the  bottom  scattering  strength  from  the 
two  very  different  bottom  scattering  kernels  (and  hence 
mechanisms)  are  similar.  Since  the  depression  angle  of  the  sonar  in 
this  study  is  3°,  the  grazing  angle  in  the  study  site  is  about  2.9°. 
Figure  47  suggests  that  there  should  be  less  than  a  5  dB  difference 
between  the  two  bottom  scattering  strength  kernels.  However,  even 
small  variations  in  the  temperature  of  the  water  mass  for  Run  18 
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could  further  decrease  the  bottom  grazing  angle  and  reduce  the 
difference  the  between  bottom  scattering  strength  kernels  (Figure 
47).  Therefore,  the  similarity  of  modeled  RL  curves  using  the  two 
different  bottom  scattering  kernels  is  not  surprising.  The  very  good 
curve  fit  of  the  two  modeled  RLs  to  the  measured  RL  is  also 
expected.  However,  the  near  coincidence  of  the  two  modeled  RL 
curves  was  not  expected  and  is  likely  due  to  the  water  temperature 
(downward  refracting)  and  hence  very  shallow  grazing  angles. 

Because  of  the  low  depression  angle  of  the  sonar,  and  hence 
low  grazing  angle  of  the  acoustic  energy,  at  least  two  bottom 
scattering  kernels  can  describe  the  bottom  scattering  strength. 

Each  of  the  kernels  represents  different  scattering  mechanisms 
governed  by  different  physics.  But  only  one  mechanism  is  correct. 
Wong  and  Chesterman  (1968)  concluded  that  for  grazing  angles 
between  0.5°  and  6°,  bottom  scattering  strength  was  consistent  for 
powers  of  sin  (0)  ranging  from  0  to  1.5  .  The  National  Defense 

Research  Council  (1946)  and  McKinney  and  Anderson  (1964)  arrived 
at  the  same  conclusion.  McCammon  (1993)  suggests  that  this  may  be 
due  to  volume  scattering  at  such  low  grazing  angles.  However, 
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constancy  of  bottom  scattering  strength  at  low  grazing  angles 
seems  more  likely  a  remnant  of  the  mathematics  (more  than  two 
free  parameters)  used  to  quantify  bottom  scattering  (10  log  p  and 
the  power  of  sin  (6))  and  the  TL  characteristic  of  the  bottom 
(Figures  41,  42,  and  43). 

There  are  then  four  possible  bottom  scattering  strength  kernel 
variations  which  fit  the  GSM  RL  estimates  to  the  measured  RL  for 
beam  1  of  ping  2  (Table  4).  Based  again  on  geology  (Figures  6  and  7), 
kernels  1  and  2  of  Table  4  are  likely  candidates  for  the  one  correct 
bottom  scattering  strength  kernel  since  beam  1  initially  propagates 
over  sand  and  then  propagates  over  silt-clay.  Considering  that 
acoustic  energy  penetrates  the  silt-clay  sediment  and  is  refracted 
in  the  sub-bottom  and  not  simply  reflected  from  the  water/sediment 
interface  (Figure  19  and  20),  it  is  not  logical  to  presume  a  simple 
Lambert’s  law  scattering  mechanism  throughout  the  propagation 
path  since  this  law  is  based  on  scattering  at  the  water/sediment 
interface.  A  volume  scattering  mechanism  (sin  (0))  is  more 
applicable  to  a  silt-clay  bottom  because  of  increased  penetration  of 
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the  acoustic  energy.  Volume  scattering  from  silt-clays  is  supported 
by  the  works  of  Jackson  and  Briggs  (1992)  and  Mourad  and  Jackson 
(1993).  This  supports  the  conclusion  that  in  the  first  7.5  seconds  of 
reverberation  for  ping  2,  beam  1,  propagation  is  over  the  sand 
bottom.  Beyond  7.5  seconds  reverberation  propagation  is  over  the 
silt-clay  bottom  type.  Hence,  kernel  number  1  of  Table  4  is  then  the 
correct  combination  of  bottom  scattering  strength  kernels  to 
quantify  and  describe  the  bottom  scattering  along  beam  1  of  ping  2. 

Thus,  even  for  a  single  beam  of  a  ping,  at  least  two 
geoacoustic  models  with  different  bottom  scattering  strength 
kernels  are  required  to  model  the  ocean  bottom  in  this  local  of  Area 
Foxtrot. 

Interestingly,  Cole  and  Podeszwa  (1974)  published  bottom 
scattering  strength  measurements  in  Area  Foxtrot  over  a  reportedly 
sand  bottom  for  grazing  angles  near  15  °.  They  showed  the  bottom 
scattering  strength  to  be  -29  +  10  log  sin  (0).  This  measurement  is 
comparable  to  the  bottom  scattering  strength  of  -35  +  10  log  sin  (0) 
for  the  silt-clay  bottom  as  determined  in  this  study. 
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Hence,  without  some  prior  detailed  knowledge  of  the  geology  of 
Area  Foxtrot,  one  could  have  easily  used  a  silt-clay  geoacoustic 
model  and  applied  a  bottom  scattering  strength  using  sin  (0)  to  the 
entire  RL  time  series  for  beam  1  of  ping  2.  This  would  have  been 
wrong  in  light  of  the  local  geology. 
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Table  4.  Four  possible  bottom  scattering  strength  kernel  scenarios 
that  fit  GSM/MPE  RL  to  measured  RL. 


Sb  kernel  # 

Time  0  -  7.5  (sec) 

Time  7.5  -  20  (sec) 

1 

-33  +10  log  sin2  (0) 

-35  +10  log  sin  (0) 

2 

-33  +10  log  sin2  (0) 

-17  +10  log  sin2  (0) 

3 

-35  +10  log  sin  (0) 

-35  +10  log  sin  (0) 

4 

-17  +10  log  sin2  (0) 

-17  +10  log  sin2  (0) 
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D.  CURVE  FITTING  OF  BOTTOM  SCATTERING  STRENGTH 

Figure  48  shows  the  two  bottom  scattering  angular 
dependencies  of  sin  (0)  and  sin2  (0)  but  neglects  the  10  log  |i  and  the 

TL  dependence  of  bottom  scattering  strength.  Volume  scattered 
energy  is  affected  by  the  TL  in  the  sediment  layer  since  it 
propagates  within  the  sediment.  Normally,  one  would  expect  higher 
RL  when  the  angular  dependence  is  sin  (0)  (Figure  48).  However,  the 
TL  within  the  sediment  can  reduce  RL  to  levels  comparable,  or  even 
less  than,  when  the  angular  dependence  is  sin2  (0).  Figures  46  and 

47  include  the  10  log  |i  variable  (10  log  |i  =  -35  and  -17)  factored 
into  the  bottom  scattering  strength  angular  dependencies  of  sin  (0) 
and  sin2  (0),  respectively.  These  figures  show  that  at  very  shallow 

grazing  angles  the  two  bottom  scattering  strength  kernels  yield 
bottom  scattering  strengths  within  only  a  few  decibels  of  each 
other.  Figures  49,  50,  51,  and  52  also  illustrate  the  similarity  of 
bottom  scattering  strengths  for  the  two  angular  dependencies  when 
different  values  of  10  log  p  are  substituted  into  the  scattering 
strength  equation.  Notice  that  as  the  difference  between  the 
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10  log  p,  values  approaches  zero,  the  bottom  scattering  strength 
curves  increase  in  similarity  over  a  wider  range  of  grazing  angles. 
Figures  51  and  52  show  that  for  the  bottom  scattering  strength 
kernels  of  -25  +  10  log  sin  (0)  and  -17  +  10  log  sin2  (e),  the  two 

bottom  scattering  strength  curves  are  within  about  2.5  dB  of  each 
other  for  grazing  angles  of  5°  to  16°.  The  curves  have  a  crossover 
point  at  about  9°.  Additionally,  as  the  values  of  10  log  p  in  the  two 

bottom  scattering  strength  kernels  approach  equality,  the  crossover 
point  for  the  two  curves  shifts  toward  higher  grazing  angles  and  the 
two  bottom  scattering  strength  curves  show  more  similarity  over  a 
wider  range  of  angles.  This  implies  that  selecting  the  correct 
angular  dependence  when  curve  fitting  measured  bottom  scattering 
strength  data  is  extremely  difficult  without  detailed  geoacoustic 
models  to  test  the  theoretical  results.  With  two  free  parameters, 

10  log  p  and  sinx  (6),  there  are  physically  meaningless  ways  to  curve 

fit  bottom  scattering  strength  data  with  “false”  accuracy. 

Wrong  assumptions  or  assumptions  based  on  little  geologic  and 
geoacoustic  data  can  easily  lead  one  to  believe  the  correct  bottom 
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scattering  strength  kernel  has  been  selected  (when  it  has  not  been). 
For  example,  if  TL  based  on  an  inaccurate  geoacoustic  model  were 
underestimated  by  about  2  dB  for  grazing  angles  less  than  9°  and  TL 
were  overestimated  by  the  same  amount  for  grazing  angles  greater 
than  9°,  the  -17  +  10  log  sin2  (e)  curve  in  Figure  52  could  be 

indistinguishable  from  the  -25  +  10  log  sin  (0)  curve.  This  could 

easily  occur  in  a  range-dependent  shallow  water  environment  and 
re-iterates  the  necessity  of  developing  accurate  geoacoustic  models 
based  on  detailed  geologic  surveys  when  studying  bottom  scattering 
strength.  The  wide-range  of  10  log  p  values  applied  to  sediment 
types  thought  to  be  the  same  (see  McCammon,  1993)  is  probably  due 
to  one  or  more  of  the  factors  discussed  above. 

Bottom  scattering  strength  and  TL  are  both  functions  of 
geology  and  sediment  type.  The  variability  of  sediment  type  (and 
thus  TL)  in  shallow  water  is  manifested  in  variability  of  bottom 
scattering  strength  kernels.  Hard  sand  bottoms  have  bottom 
scattering  kernels  that  obey  Lambert’s  law  while  silt-clay  bottoms 
show  more  volume  scattering  and  do  not  obey  Lambert’s  law. 
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VI.  CONCLUSIONS 


The  current  geoacoustic  model  for  Area  Foxtrot  was  found  to 
be  inaccurate  because  the  sediment  compressional  wave  speed  and 
density  and  attenuation  values  did  not  conform  to  the  geoacoustic 
modeling  theories  given  by  Hamilton  and  Bachman  (1982).  The 
surface  sediments  in  Area  Foxtrot  are  also  much  more  spatially 
variable  than  originally  thought.  New  Hamilton-type  geoacoustic 
models  were  developed  for  the  area.  To  calculate  RL,  the  FEPE 
propagation  model  was  used  to  augment  the  GSM  propagation  model 
(currently  used  by  NUWC  for  TL  estimation  and  RL  modeling)  because 
the  bottom  loss  model  in  GSM  did  not  estimate  TL  accurately  in 
shallow  water.  The  comparison  between  modeled  RL  and  measured 
RL  (from  a  1991  ASW  exercise  conducted  by  NUWC)  enabled  bottom 
scattering  strength  kernels  to  be  developed  for  Area  Foxtrot.  The 
conclusions  reached  in  this  study  are  as  follows: 

1 )  Based  on  geologic  data,  it  is  readily  apparant  that  Area  Foxtrot 
is  not  a  spatially  homogeneous,  acoustically  benign  area.  The 
surface  sediment  was  found  to  range  from  course  sand  and  shell 
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fragments  to  silt-clay  sediments.  Hamilton-type  geoacoustic 
models  developed  in  this  study  show  that  compressional  wave 
speeds  in  Area  Foxtrot  vary  from  over  1800  m/s  to  less  than  1600 
m/s  over  ranges  of  less  than  7  km.  No  one  geoacoustic  model  can 
describe  all  of  Area  Foxtrot. 

2)  The  current  geoacoustic  model  for  Area  Foxtrot  is  not 
consistent  with  either  observed  geological  data  or  with  Hamilton- 
type  geoacoustic  modeling  theories.  The  current  “sand”  bottom 
model  for  Area  Foxtrot  employs  a  silt-clay  compressional  sound 
speed  but  uses  a  density  and  attenuation  value  consistent  with  a 
sand  sediment. 

3)  The  GSM  model  does  not  have  a  sub-routine  to  incorporate  sub¬ 
bottom  geoacoustic  properties.  A  second,  external  TL  model  must  be 
used  to  accurately  estimate  TL  in  shallow  water  environments 
where  bottom  interaction  is  of  major  concern.  NUWC  currently  uses 
bottom  loss  curves  generated  by  the  SAFARI  model  to  augment  GSM. 
NFS  uses  FEPE  TL  estimates  to  address  this  problem  because  FEPE 
has  repeatedly  demonstrated  an  ability  to  accurately  model  sub- 
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bottom  refracted  energy.  Either  way,  GSM  can  not  model  bottom 
interacting  energy  accurately  without  using  another  TL  model. 

4)  Significant  differences  were  noted  in  the  TL  estimates 
between  GSM  and  SAFARI  or  FEPE  for  both  the  sand  and  silt-clay 
geoacoustic  models.  These  differences  in  turn  significantly 
influenced  the  shape  of  the  modeled  RL  curve  demonstrating  the 
sensitivity  of  RL  to  TL. 

5)  Because  the  bottom  surface  sediments  in  Area  Foxtrot  vary 
from  course  sand  to  silt-clay,  the  bottom  scattering  strength 
kernels  vary  also.  The  sand  bottom  was  well-modeled  using  a 
bottom  scattering  strength  of  -33  +  10  log  sin2  (0)  and  is  consistent 
with  a  Lambert’s  law  relationship.  The  silt-clay  bottom 
demonstrated  a  bottom  scattering  strength  of  -35  +  10  log  sin  (0) 

and  is  consistent  with  the  volume  scattering  mechanism  proposed  by 
Jackson  and  Briggs  (1992)  and  Mourad  and  Jackson  (1993). 

6)  When  the  GSM  TL  model  was  corrected  for  bottom  loss,  i.e., 
incorporating  bottom  refraction  as  modeled  by  FEPE,  the  following 
results  were  obtained: 
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a.  Excellent  agreement  was  achieved  between  modeled 
RL  and  measured  RL  for  beams  steered  in  the  direction  of  the  hard 
sand  bottom  using  Lambert’s  law  (sin2  (e)  dependence)  to 

characterize  the  bottom  scattering  strength. 

b.  Excellent  agreement  was  achieved  between  modeled  RL  and 
measured  RL  for  beams  steered  in  the  direction  of  the  silt-clay 
bottom  using  a  more  “lossy”  (sin  (0)  dependence)  characterization  of 

the  bottom  scattering  strength.  This  is  consistent  with  volume 
scattering  as  described  by  Jackson  and  Briggs  (1992)  and  Mourad  and 
Jackson  (1993). 

c.  More  than  one  bottom  scattering  strength  kernel  can  be  used 
to  fit  the  modeled  RL  for  the  same  bottom  type  to  the  measured  RL 
for  very  shallow  grazing  angles  (less  than  3°).  As  the  difference  in 
the  10  log  p  terms  approaches  zero  in  two  different  bottom 

scattering  strength  kernels  (with  angular  dependencies  of  sin  (0)  and 
sin2  (0)),  the  resulting  bottom  scattering  strength  curves  increase  in 
similarity  over  a  wider  range  of  grazing  angles. 
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7)  The  performance  of  the  AN/SQS-53C  (CW  500  pulse)  in 
shallow  water  is  strongly  dependent  on  the  local  geology  of  the 
bottom  in  negative  (downward)  refracting  water  columns.  The  RL 
will  fall  off  more  quickly  when  propagation  is  over  a  silt-clay 
bottom  than  over  a  hard  sand  bottom. 

8)  The  modeling  techniques  developed  in  this  study  can  be  used  to 
bound  RL  in  active  sonar  performance  in  strategic  shallow  water 
areas  where  little  is  known  about  the  bottom  sediment  type  and  its 
spatial  distribution.  As  shown  in  the  bottom  panel  of  Figure  41,  the 
RL  values  can  be  modeled  for  hard  sand  and  silt-clay  sediments 
based  on  Hamilton  modeling  principles  and  by  using  the  appropriate 
angular  dependence  of  bottom  scattering  strength  (sin2  (e)  for  sand, 

sin  (0)  for  silt-clay).  Measured  RL  should  then  fall  between  these 
two  extremes  for  locations  where  the  bottom  is  sediment  covered. 

A  method  can  easily  be  developed  in  which  measured  RL  from  a  few 
pings  emitted  in  the  area  of  interest  could  be  used  to  deduce  the 
azimuthal  dependence  of  RL.  By  knowing  the  azimuthal  dependence 
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of  RL  in  the  local  area,  the  bottom  sediment  type  can  be  deduced  and 
more  accurate  detection  ranges  can  be  determined. 

9)  Bottom  scattering  strength  is  a  function  of  grazing  angle, 
Mackenzie’s  constant,  and  TL.  The  TL  must  be  accurately  modeled  to 
deduce  the  bottom  scattering  strength  especially  at  low  grazing 
angles  for  which  volume  reverberation  within  the  sediment  layer  is 
important.  Incorrect  bottom  scattering  strength  kernels  can  easily 
be  deduced  from  curve  fitting  to  measured  data.  This  occurs  because 
various  combinations  of  Mackenzie’s  constant  (10  log  p),  first  or 
second  power  angular  dependence  (sin  (0)  or  sin2  (0)),  and  TL  trends 
(due  to  propagation  over  different  bottom  sediments)  can  produce 
very  similar  RL  curves. 
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VII.  RECOMMENDATIONS 


The  1992  Navy  and  Marine  Corps  White  Paper  “...From  the  Sea” 
authored  by  the  Secretary  of  the  Navy,  the  Chief  of  Naval  Operations, 
and  the  Commandant  of  the  Marine  Corps  (O’Keefe  et  al.,  1992) 
states  “the  shift  in  strategic  landscape  means  Naval  Forces  will 
concentrate  on  littoral  warfare  and  maneuver  from  the  sea.” 
Maneuvering  from  the  sea  “provides  a  potent  war-fighting  tool  to 
the  Joint  Task  Force  Commander-a  tool  that  is  literally  the  key  to 
success  in  many  likely  contingency  scenarios.” 

The  following  recommendations  derived  from  the  results  of 
this  thesis  are  provided  in  support  of  the  War-fighting  concepts 
envisioned  in  “...From  the  Sea.”  These  recommendations  are  designed 
to  improve  the  accuracy  of  performance  prediction  of  current  and 
emerging  active  sonar  systems  and  put  forth  ideas  for  the  tactical 
applications  of  this  thesis. 

A.  RECOMMENDATIONS  REGARDING  AREA  FOXTROT 

1)  If  NUWC  or  the  Navy  plan  to  use  Area  Foxtrot  as  a  sonar  test  bed 
for  active  sonar  systems,  then  it  is  strongly  recommended  to 
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conduct  a  detailed  geologic  study  of  at  least  two  15  km  by  15  km 

sections  of  Area  Foxtrot.  One  section  should  have  a  predominantly 

hard  sand  bottom  while  the  other  has  a  silt-clay  bottom.  The 

geologic  study  must  include  surface  sediment  samples,  sediment 

cores,  and  bathymetry.  The  resolution  of  the  geologic  survey  will 

determine  the  accuracy  of  any  resulting  geoacoustic  models. 

2)  Construct  Hamilton-type  geoacoustic  models  for  the  sand  and 

silt-clay  bottom  areas  and  use  these  models  to  improve  the  accuracy 

of  performance  prediction  of  Navy  active  sonars  in  Area  Foxtrot. 

B.  RECOMMENDATION  REGARDING  BOTTOM  SCATTERING 
STRENGTH 

1)  Use  a  bottom  scattering  strength  kernel  with  angular  dependence 
of  sin2  (e)  for  hard  sand  sediment  type  bottoms  and  sin  (0)  for  silt- 

clay  sediment  bottoms.  Mackenzie’s  constant  (10  log  p)  Is  variable 
in  both  sand  and  silt-clay  sediments  and  should  be  determined  by 
curve  fitting  modeled  RL  to  measured  RL  while  applying  the 
appropriate  angular  dependence  for  the  sediment  type. 
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C.  RECOMMENDATION  REGARDING  GSM 


1)  Integrate  FEPE  into  the  GSM  model  for  estimating  TL  Since  FEPE 

can  take  several  hours  to  run,  it  is  not  recommended  for  tactical 

operational  use.  For  tactical  applications,  a  database  of  FEPE  TL 

curves  can  be  calculated  in  the  laboratory  now  for  shallow  water 

environments  of  interest  (or  potential  conflict)  based  on  Hamilton 

geoacoustic  models  developed  for  these  areas.  The  TL  curves  can  be 

maintained  in  a  library  within  GSM  and  accessed  by  latitude, 

longitude  and  ship’s  heading.  This  would  permit  a  shipboard  METOC 

(Meteorology  /  Oceanography)  officer  to  get  a  first-cut  estimate  of 

TL  and  RL  in  the  operating  area  based  on  full  Hamilton  geoacoustic 

models.  Recall  that  FEPE  is  range-dependent  so  TL  curves  (and  hence 

RL  curves)  can  be  calculated  as  range-dependent.  Development  of 

geoacoustic  models  by  the  METOC  officer  is  not  recommended. 

D.  RECOMMENDATION  REGARDING  THE  TACTICAL 
APPLICATION  OF  THIS  THESIS 

1)  Develop  the  method  suggested  in  this  study  where  a  few 

(nominally  12)  on-site  RL  curves  can  be  used  to  deduce  the  gross 

bottom  sediment  type  in  an  operating  area  based  on  the  measured  RL 
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curve  shapes  (see  the  flow  chart  in  Figure  53).  Information  on 
sediment  type  can  be  input  to  GSM  to  make  the  laboratory-developed 
geoacoustic  model  more  accurate.  This  would  certainly  improve  the 
accuracy  of  GSM  TL  and  RL  estimates  and  in  turn  provide  more 
accurate  range  predictions  in  tactical  environments  where  little  is 
known  about  the  local  geoacoustic  environment. 
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VIII.  APPENDIX.  FIGURES 
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SOUND  VEIOCITY,  m/i 


Figure  1.  Porosity  versu 
samples  from  the  contini 
Bachman,  1982). 


MEAN  GRAIN  SIZE 


Figure  2.  Mean  grain  size  versus  compressional  wave  speed 
in  samples  from  the  continental  terrace  (from  Hamilton 
Bachman,  1982). 


104 


Attenuation  of  Compressional  Waves  (kp) 
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Figure  3.  Compressional  wave  attenuation  versus  depth  for  sand 

and  silt-clay  sediments  on  the  continental  terrace  (from  Hamilton,  1980). 
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Sound  Speed  (m/s) 


Figure  4.  Compressional  wave  speed  versus  depth  from  the 
water  surface  to  the  bottom  of  the  sediment  layer  -  the  current 
geoacoustic  model. 
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Figure  5.  Locations  of  sediment  samples  taken  in  Area  Foxtrot. 
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Figure  6.  Contour  plot  of  surface  sediment  mean  grain  size 
(phi  units)  in  Area  Foxtrot. 
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Figure  7.  Contour  plot  of  compressional  wave  speed  calculated  from 
Hamilton  and  Bachman’s  (1982)  regression  equation.  No  correction  is 
made  for  geology  or  standard  error  in  the  calculation. 


Approximate  location  and  Approximate  location  and  track 

track  of  the  surface  ship.  of  the  target. 


Runs  11-18  Runs  11-18 


Figure  8.  Same  as  Figure  7  with  the  approximate  locations  and  tracks 
of  the  surface  ship  and  target  submarine  for  the  3-13  September  1991 
shallow  water  ASW  test  overlain. 


no 


Latitude  -  North 


Longitude  -  West 

Figure  9.  Surface  sediment  parameters  surrounding  the  Run  18  area. 
Black  diamonds  represent  the  location  of  a  sediment  sample. 

Legend,  for  each  location;  mean  grain  size  in  phi  units 

grain  constituents  in  %  as  gravel/sand/silt/clay 

compressional  wave  speed  in  m/s 

lithology  code:  A  -  Green  silty  sand  and  shells 

B  -  Coarse,  brown  well-sorted  sand 

C  -  Green  silty  sand  and  shell  debris 

D  -  Green  silt  and  shell  fragments 

E  -  Brown  to  greyish  medium  to  coarse  sand 

F  -  Coarse  brown  sand 

G  -  Brownish-green  very  sandy  clay 
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71.6 


0 


Sound  speed  (m/s) 


Figure  10.  Compressional  wave  speed  versus  depth  from 
the  water  surface  to  the  bottom  of  the  sediment  layer  - 
NPS  geoacoustic  model. 
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Figure  11.  To  calculate  bottom  reverberation  TL  was  calculated 
from  the  source  to  the  ocean  bottom  as  the  receiver  on  the  bottom 
was  moved  out  in  range. 
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Figure  12  A  plot  of  TL  from  the  GSM/MPE  model  using  the  NFS  sand  bottom  model  as  input. 

The  yellow  line  is  random  (incoherent)  TL,  the  white  line  is  phased  (coherent)  TL,  the  blue  line 
represents  spherical  spreading,  and  the  red  line  is  direct  path  TL.  The  eigenray  tolerance 
is  0.1.  The  sediment  attenuation  in  the  plot  legend  is  in  nepers/wave  . 


FigureTs^  A  plot  of  TL  from  the  GSM/MPE  model  using  the  current  bottom  model  as  input.  The 
yellow  line  is  random  (incoherent)  TL,  the  white  line  is  phased  (coherent)  i  L,  the  blue  line 
represents  spherical  spreading,  and  the  red  line  is  direct  path  TL.  The  eigenray  tolerance  is 
0.1.  The  sediment  attenuation  in  the  plot  legend  is  in  nepers/wave  . 


Figure  14.  Same  as  Figure  12  except  the  eigenray  tolerance  is  0.001. 


Paitss 


Figure  16.  Propagation  loss  using  the  FEPE  TL  model  for  the  NPS  sand  bottom  model  (red 
curve)  and  the  current  bottom  modei  (yellow  curve).  The  source  is  at  7.3  m  and  the  receiver  is 
at  6  4  m  (on  the  bottom). 


Figure  17.  Same  as  Figure  16  except  the  source  and  receiver  are  both  at 


Figure  19.  Comparison  between  the  GSM/MPE  TL  model  (the  yeiiow  curve  is  phased  TL,  the 
blue  curve  is  random  TL)  and  the  FEPE  TL  model  (red  curve)  using  the  current  geoacoustic 
model  as  input.  This  illustrates  that  GSM/MPE  underestimates  the  propagation  loss  compared  to 
the  FEPE.  GSM/MPE  uses  the  Rayleigh  reflection  coefficient  for  bottom  loss.  FEPE  includes  a 
fuii  Hamilton-type  geoacoustic  model  for  TL  estimation. 


Figure  20.  Same  as  Figure  19  except  GSM/MPE  phased  TL  is  not  included. 


Figure  21.  Comparison  between  the  GSM/MPE  TL  model  (the  yellow  curve  is  phased  TL,  the 
blue  curve  is  random  TL)  and  the  FEPE  TL  model  (red  curve)  using  the  NPS  sand  model  as  input. 
This  illustrates  that  GSM/MPE  overestimates  the  propagation  loss  compared  to  the  FEPE. 
GSM/MPE  uses  the  Rayleigh  reflection  coefficient  for  bottom  loss.  FEPE  includes  a  full 
Hamilton-type  geoacoustic  model  for  TL  estimation. 


Figure  22.  Same  as  Figure  21  except  GSM/MPE  phased  TL  is  not  included 


Figure  23.  Derived  backscattering  strength  at  100  Hz 
(from  Hanna  et  al.,  1994). 
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Figure  24.  Derived  backscattering  strength  at  800  Hz 
(from  Hanna  et  al.,  1994). 
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Figure  25.  Relative  loss  associated  with  low  grazing  angle  paths 
and  high  grazing  angle  paths  (from  Hanna  et  al.,  1994). 
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Figure  26.  Measured  RL  and  GSM/MPE  RL  with  Sb  set  to  zero  for  ping  2 
beam  1  (top  panel).  Measured  RL  minus  GSM/MPE  RL  with  a  3rd-order 
polynomial  cure  fit  (bottom  panel). 
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Figure  27.  Measured  RL  and  GSM/MPE  RL  with  Sb  set  to  zero  for  ping  2, 
beam  12  (top  panel).  Measured  RL  minus  GSM/MPE  RL  with  a  3rd-order 
polynomial  cure  fit  (bottom  panel). 
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Figure  28.  Measured  RL  and  GSM/MPE  RL  with  Su  set  to  zero  Pitg  2 
beam  24  (top  panel).  Measured  RL  minus  GSM/MPE  RL  with  a  3rd-order 
polynomial  cure  fit  (bottom  panel). 
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Figure  29a.  Measured  RL  and  GSM/MPE  RL  using  the  NFS  model  for  ping  2, 
beam  1  (top  panel),  ping  2,  beam  3  (center  panel),  and  ping  2,  beam  5 

(bottom  panel). 
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Figure  29b.  Measured  RL  and  GSM/MPE  RL  using  the  NFS  model  for  ping  2 
beam  7  (top  panel),  ping  2,  beam  9  (center  panel),  and  ping  2,  beam  11 
(bottom  panel). 
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Figure  29c.  Measured  RL  and  GSM/MPE  RL  using  the  NPS  model  for  ping  2, 
beam  13  (top  panel),  ping  2,  beam  15  (center  panel),  and  ping  2,  beam  17 
(bottom  panel). 
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dB// 1  uPa/Hz^1/2 


Figure  29d.  Measured  RL  and 
beam  19  (top  panel),  ping  2, 


GSM/MPE  RL  using  the  NPS  model  for  ping  2, 
beam  21  (center  panel),  and  ping  2,  beam  23 


(bottom  panel). 
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dB  // 1  uPa  /  Hz''1/2 


Figure  30a.  Measured  RL  and  GSM/MPE  RL  using  the  NP3  model  for  ping  30. 
beam  1  (top  panel),  ping  30,  beam  3  (center  panel),  and  ping  30,  beam  5 
(bottom  panel). 
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Figure  30c.  Measured  RL  and  GSM/MPE  RL  using  the  NPS.  model  for  ping  30, 
beam  13  (top  panel),  ping  30,  beam  15  (center  panel),  and  ping  30,  beam  17 
(bottom  panel). 
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Figure  30d.  Measured  RL  and  GSM/MPE  RL  using  the  NPS.  model  for  ping  30, 
beam  19  (top  panel),  ping  30,  beam  21  (center  panel),  and  ping  30,  beam  23 
(bottom  panel). 
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Figure  31b.  Measured  RL  and  GSM/MPE  RL  using  the  NPS  model  for  ping  60, 
beam  7  (top  panel),  ping  60,  beam  9  (center  panel),  and  ping  60,  beam  11 
(bottom  panel). 
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PING;  2  Measured  Reverberation  Levels  in  dB 


PING:  20  Measured  Reverberation  Levels  in  dB 


Figure  33.  Same  as  Figure  32  except  for  ping  20. 
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Figure  34.  Same  as  Figure  32  except  for  ping  40. 


PING:  60  Measured  Reverberation  Levels  in  dB 
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Figure  35.  Same  as  Figure  32  except  for  ping  60. 


PING:  88  Measured  Reverberation  Levels  In  dB 
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Figure  36.  Same  as  Figure  32  except  for  ping  88. 
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Figure  37.  Same  as  Figure  36  except  plotted  in  3-D. 


dB  // 1  uPa  /  HzA1/2  dB  // 1  uPa  /  Hz^l/a 


Figure  38a.  Comparison  between  the  Measured  RL  and  GSM/MPE  RL 
(without  corrected  TL)  using  the  NPS  silt-clay  model  as  input  (top  panel) 
for  ping  2,  beam  1.  Measured  RL  and  GSM/MPE  RL  (with  corrected  TL)  using 
three  different  values  of  10  log  p  (bottom  panel).  The  angular  dependence 
is  modeled  as  sin  (0). 
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Figure  38b.  Same  as  Figure  38a  except  for  ping  20,  beam  1 


dB  // 1  uPa  /  Hz^1/2  dB  // 1  uPa  /  Hz^1/2 


Figure  38c.  Same  as  Figure  38a  except  for  ping  40,  beam  1 . 
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Figure  38e.  Same  as  Figure  38a  except  for  ping  88,  beam  1 
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Figure  39a.  Comparison  between  the  Measured  RL  and  GSM/MPE  RL 
(without  corrected  TL)  using  the  NPS  silt-clay  model  as  input  (top  panel) 
for  ping  2,  beam  12.  Measured  RL  and  GSM/MPE  RL  (with  corrected  TL) 
using  three  different  values  of  10  log  p  (bottom  panel).  The  angular 
dependence  is  modeled  as  sin  (0). 
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Figure  39b.  Same  as  Figure  39a  except  for  ping  20,  beam  12. 


Measured  RL 

GSM/MPE  RL  w/o  corrected  TL 

_ , - -  I  I  I _ L _ L. 

4  6  8  10  12  14  16 

Time  (sec) 


T 


- Measured  RL 

- (Sb  =  -30  +  10  log  Sin^2  (theta)) 

- (Sb  =  -35  +  10  log  Sin^2  (theta)) 

(Sb  =  -40  +  10  log  Sin^2  (theta)) 


w/  corrected  TL 


,  ,  I  I - '  I 

2  4  6  8  10  12 

Time  (sec) 


_l _ L_ 

14  16 


Figure  39c.  Same  as  Figure  39a  except  for  ping  40,  beam  12 
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dB  // 1  uPa  /  Hz^1/2  dB  // 1  uPa  /  Hz^1/2 


Figure  39e.  Same  as  Figure  39a  except  for  ping  88,  beam  12. 
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Figure  40a.  Comparison  between  the  Measured  RL  and  GSM/MPE  RL 
(without  corrected  TL)  using  the  NPS  silt-ciay  model  as  input  (top  panel) 
for  ping  2,  beam  24.  Measured  RL  and  GSM/MPE  RL  (with  corrected  TL) 
using  three  different  values  of  10  log  p,  (bottom  panel).  The  angular 
dependence  is  modeled  as  sin  (6). 
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dB  // 1  uPa  /  Hz^1/2  dB  // 1  uPa  /  Hz'M/2 


Figure  40b.  Same  as  Figure  40a  except  for  ping  20,  beam  24. 
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Figure  40d.  Same  as  Figure  40a  except  for  ping  60,  beam  24 
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Figure  40e.  Same  as  Figure  40a  except  for  ping  88,  beam  24. 
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Fioi.rp  41  Comoarison  between  the  measured  RL  and  GSM/MPE  RL  (without  corrected  TL)  for 
So  2  beamT  us  ng  the  sand  model  (red  curve)  and  NPS  silt-clay  model  (blue  curve)  as 
S  (iop  pTnel').  comparison  bstoean  the  measured  RL  and  GSM/MPE  RL  JU 

us^nq  the  NPS  sand  model  with  bottom  scattering  strength  angular  dependence  of  ^ 

curve)  and  NPS  silt-clay  model  with  bottom  scattering  strength  angular  ^ 

(blue  curve)  (bottom  panel).  Notice  that  RL  for  the  two  different  bottom 

bottom  scattering  angular  dependencies  (but  similar  10  log  p)  fit  the  measured  RL  out  to  about 

7.5  seconds  after  which  the  curves  diverge. 


dB  // 1  uPa  /  Hz^1/2  dB  // 1  uPa  /  Hz^1/2 


Figure  42.  Same  as  Figure  41  except  that  the  bottom  scattering  strength 
kernels  are  interchanged  between  the  two  bottom  types  to  essentially 
separate  the  TL  trend  in  the  bottom  type  models  from  the  angular 
dependencies.  Notice  now  that  the  two  modeled  RLs  show  no  similarity. 
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Figure  43.  Measured  RL  (green  curve)  and  modeled  RL  using  the  NPS  silt- 
clay  bottom  as  input  with  two  different  bottom  scattering  strength 
kernels;  -17  +  10  log  sin2  (e)  (red  curve)  and  -35  +  sin  (0)  (blue  curve). 
Note  that  the  two  modeled  RL  curves  are  nearly  coincident  throughout  the 
entire  time  series. 
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Figure  44.  Expanded  view  of  Figure  43. 
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Figure  45.  Same  as  Figure  43  except  that  modeled  RL  using  the  NFS  sand 
bottom  with  bottom  scattering  kernel  of  -33  +  10  log  sin2  (6)  is  included 
(light  blue  curve).  This  illustrates  the  similarity  of  the  three  modeled 
RLs  (with  different  bottom  scattering  strength  kernels)  out  to  about  7.5 
seconds. 


Bottom  Scattering  Strength  (dB) 


Figure  46.  Bottom  scattering  strength  as  a  function  of  grazing  angle  for 
the  bottom  scattering  strength  kernels  -35  +  10  log  sin  (0)  (red  curve)  and 
-17  +  10  log  sin2  (0)  (green  curve).  The  difference  between  the  curves  is 
small  at  very  low  grazing  angles. 


Grazing  Anglo 


Figure  48.  Bottom  scattering  strength  as  a  function  of  grazing  angle  to 
illustrate  the  scattering  angular  dependencies,  sin  (6)  and  sin2  (0).  Note 
that  the  10  log  p  term  is  neglected  in  the  calculation  of  bottom  scattering 
strength  (from  McCammon,  1993). 
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Figure  49.  Bottom  scattering  strength  as  a  function  of  grazing  angle  for 
the  bottom  scattering  strength  kernels  -30  +  10  log  sin  (9)  (red  curve)  and 
-17  +  10  log  sin2  (6)  (green  curve).  Notice  that  the  curves  show  more 
similarity  as  the  difference  between  the  two  10  log  |a.  terms  is  reduced 
(compare  to  Figure  46). 
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Figure  51.  Bottom  scattering  strength  as  a  function  of  grazing  angle  for 
the  bottom  scattering  strength  kernels  -25  +  10  log  sin  (6)  (red  curve)  and 
-17  +  10  log  sin2  (0)  (green  curve).  Notice  that  as  the  difference  between 
the  two  10  log  |a.  terms  is  reduced  (from  that  in  Figures  46  and  49)  the 
bottom  scattering  strength  curves  increase  in  similarity  and  remain 
similar  over  a  wider  range  of  grazing  angles. 
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